AO-A043  547 
UNCLASSIFIED 


GENERAL  RESEARCH  CORP  SANTa  BARBARA  CALIF  F/6  l7/o 

2®SC0E  manual*  volume  17,  HIGH-ALTITUOE  debris-energy  DEPOS— ETC(U) 
SEP  75  0 A HAMLIN,  J Y WANG,  M R SCHOONOVER  DNA001-74-C— 0182 

DNA-3964F-17  ML 


1 or  2 

^04354? 

■ 

1 

! 

1 

t 

1 

i 

1 — ! 

t 

1 

1 

1 

) 

j 

I 

1 

I 

1 

1 

1 

1 

— 1 

_ 

_ 

V 


ir:) 

CO 


' 1 
/ 


DNA3964F-17 


THE  ROSCOE  MANUAL 


Volume  17-High  Altitude  Debris-Energy  Deposition 


^ Science  Applications,  Inc. 

<3:  P.  o.  Box  2351 

La  Jolla,  California  92038 


22  September  1975 


Final  Report  for  Period  1 March  1974—31  January  1975 
CONTRACT  No.  DNA  001 -74-C-01 82 


APPROVED  FOR  PUBLIC  RELEASE, 
DISTRIBUTION  UNLIMITED. 


n : 


c 


THIS  WORK  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODES  B322074464  S99QAXHC06428  AND 
B322075464  S99QAXHC06432  H2590D 


O- 

coprepared  for 


Director 

UJ 

- DEFENSE  NUCLEAR  AGENCY 

p * Washington,  D.  C.  20305 


UNCLASSIFIED 


security  C*.ASSlFlCATtON  OF  TmIS  PAGE  '»>!<•/»  Fnffr.  J 


, , ! REPORT  DOCUMENTATION  PAGE 

RKAD  INSTRUCTIONS 
HKKORK  COMPI-KTINr.  KORM 

“i  ^ ^ > GOVT  ACCESSION  NO 

DNA  3964F-17  ) 

3 RFCIPiEn'^'S  Ca-'-a^OG  number 

THE  ROSCOE  MANUAL  • ^ ( / 

Volume  17— High-Altitude  Debris-Energy  ^ 

Deposition 

■5  -F^Pf ’5F'"Rn*>'«,l_iJ--_^10D  COVERED 

Final /Repwt, for  i%l^d 
1 Mar  74 — 31  Jan  75,,  )' 

^ i.iRiJ  - Rfc  B E R 

SA1-75-609-LJ,  Vol.  5 

7 AjTmOR^s’  / ^ 

Daniel  A.  Hamlin.  Jon  Y.^AVang, 

Melvin  R . Schoonover  . i i 

John  I..  Valerio  ' ' 

8 contract  or  GRAN‘S  NUMBER  * 

DNA  001-74-C-0182 

^ T E RFORMING  ORGANI  Z ATION  NAME  AND  ADDRESS 

Science  Applications,  Inc. 

P.O.  Box  2351 

La  Jolla,  California  92038 

10  program  Element  PROJECT  ■'aSk 
AREA  ft  WORK  UNi"'  numbers 

NWED  Subtasks 
S99QAXHC064-28/32 

1’  TON'‘t?OLLlNG  OFFICE  NAME  AND  ADDRESS 

Director 

Defense  Nuclear  Agency 
Washington,  D.C.  20305 

12  report  DAT£ 

22  September  1975  ^ 

’ njmbFR  OF  pages  ' 

126  ! T-  ?'  r'  ' - • 

'4  MON'TQRING  agency  N AME  ft  ADDRFSSOf  T»f/cr.-nf  trom  C-mlroU ink!  O/fi.  ■ 

UNCLASSIFIED 

ISrt  DECLASSIFICATION  DOWNGRADING 

schedule 

1 Ub  DtS"^  Rl  OUTlON  ST  A*'"  EMEN  T •.,(  r/u-.  Kcpmli  I 

Approved  for  public  release:  distribution  unlimited. 


'ilSTRiBuTIOM  S‘’’ATFmFNT  'oI  'he  f enfere./  m Hl'u-k  ; / .//fh'renf  from  A’epo 


IB  supplementary  notes 


This  work  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMSS 
Codes  B322074464  S99QAXHC06428  and  B322075464  S99QAXHC06432  H2590D, 


/ idrntifv  hi  M<>.  k rmmher 

Energy  Deposition  by  Debris  Particles 
(Loss-Cone,  Ion-Leak,  and  Charge 
Exchanged 


I fi  f'l  M life  ni;rrh«r 


K F Y -VORDS  ■ ' ' -rr-Mi.,.-  f.--.  .T-  e /f  . 

ROSCOE 

Radar  Blackout  Code 
Energy  Deposition  by  X-Rays 
Energy  Depositions  by 
Ultrayi olet  Radiation 

AB'.*RA  t ‘ ..nfm’i’  ■ >n  rf  \ rr  <i  Ir  il  ru-<  t- - t 

Models  of  the  high-altitude  debris-energv  partition  and  deposition  have 
been  adopted  for  use  in  ROSCOE.  The  models  for  the  debris-energy 
partition  and  heavy-particle  source  spectra  incorporate  the  work  of 
Crevier  and  Kilb  for  the  loss-cone  and  ion-leak  particles.  The  determi- 
nation of  the  total  (and  spectrum  of  the)  UV  portion  of  the  kinetic  yield 
remaining  after  that  assigned  to  all  the  heavy -particle  motion  is  based 
on  the  work  of  Fajen  and  Sappenfield.  Representative  points  in  the  spatial 
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18.  SUPPLEMENTARY  NOTES  (ContinuecH 

EDITORS  NOTE 

\'’olumes  13  to  17  were  originally  published  by  SAI  to  describe  the  atmo- 
spheric. geomagnetic,  and  high-altitude  energy  deposition  and  neutral  heave 
models  for  ROSCOE . This  whole  section  of  code,  when  associated  with  an 
appropriate  DRIVER  subroutine,  operated  as  a package  that  ran  independently 
of  the  rest  of  the  ROSCOE  structure.  Provision  was  also  made,  within  this 
high-altitude  package,  for  two  completely  independent  descriptions  of  atmo- 
spheric heave,  each  with  its  own  description  of  atmospheric  chemistry. 

When  GRC  incorporated  this  section  of  code  within  the  ROSCOE  frame- 
work, some  modifications  were  necessary,  which  means  that  some  of  the 
descriptions  in  Volumes  13  to  17  are  inappropriate  to  ROSCOE  as  it  now 
exists.  In  particular,  the  NRL  heave  routines  (deck  NRLHYD)  and 
associated  chemistry  (deck  NRLCHM'  are  not  presently  used  in  ROSCOE. 
Three  other  subroutines  are  different ; subroutines  .AT MOSU  . EIF.  and 
XTCOEF  correspond  to  the  ROSCOE  subroutines  ATMOS,  EXPINT,  and 
WDXP  respectively.  With  these  e.xceptions.  the  subroutines  described  in 
Volumes  13  to  17  correspond  exactly  to  those  currently  in  ROSCOE. 

20.  ABSTRACT  (Continued) 

distribution  function  specified  for  the  loss-cone  and  ion-leak  particles  are 
used  as  effective  source  points  through  which  magnetic  field  lines  are  traced 
in  the  downward  direction  through  the  intercepted  grid  cells  in  the  geocentric 
quadrupole  coordinate  system.  Conventional  heavy-particle  range-energy 
theory  is  employed  along  these  paths,  without  regard  to  spiralling  effects,  to 
deposit  the  energy  of  these  heavy  particles  in  the  traversed  cells.  Charge- 
exchange  particles  are  deposited  without  regard  to  the  magnetic  field,  the 
total  heaw -particle  energies  deposited  by  inelastic  and  by  elastic  collisions  I 
in  each  cell  are  partitioned  into  various  modes  and  species  which  are  ulti-  ! 
mately  made  compatible  with  the  late-time  grid  chemistry.  The  UV  energy  j 
groups,  as  well  as  the  x-ray  energy,  are  deposited  by  tracing  ray  paths  | 

(and  constructing  line  integrals  of  relevant  species)  from  the  event  point  to 
each  cell  in  the  high-altitude  grid  treated  as  a target  cell,.-  The  portion  of  all  i 
the  deposited  energy  available  for  hydro  is  computed  as  well  as  the  ' 

instantaneously-produced  species  to  be  used  as  starting  conditions  for  the 
subsequent  hydrodynamics  and  chemistry.  Herein  are  presented  details  of 
the  models  for  tlie  debris-energy  partition  and  deposition,  including  detailed 
flow  charts  for  the  more  complicated  routines.  j 
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1.  CONTROL  ROUTINE  FOR  DEPOSITION  OF  THE  ENERGY  IN 
THE  HIGH-ALTITUDE  GRID  (SUBROUTINE  PROMPG) 


Subroutine  PROMPG,  called  by  DRIVER,  controls  the  deposi- 
tion of  the  ener^  in  the  high-altitude  grid  (HAG).  The  overall  procedure 
is  complicated,  as  Subroutine  PROMPG  directly  calls  12  subroutines 
which  in  turn  call  an  additional  26  subroutines.  The  tasks  performed  by 
Subroutine  PROMPG  include  the  following: 

a.  Tests  whether  or  not  the  event  point  is  above  the  bottom  of 
the  grid.  If  it  is,  the  program  proceeds  to  the  next  task;  if  not,  the 
program  returns  to  DRIVER  after  printing  an  appropriate  message.  It 
may  be  desirable  to  deposit  at  least  the  Group->  and  possibly  the 
Group-X  energy  for  certain  events  below  the  bottom  of  the  grid,  but  this 
feature  has  not  been  implemented  yet. 

b.  Calls  Subroutine  BOUNDY  to  set  up  the  array  GR(NCEL) 
for  the  geocentric  radius  to  the  top  boundary  of  each  cell  NCEL  in  the 
HAG.  Subroutine  BOUNDY  is  different  in  SAIHYD  than  in  NRLHTD  be- 
cause the  definitions  of  a cell  boundary  are  different  in  the  two  versions 
of  the  hydrodynamics. 

c.  Sets  up  the  array  RHOKK(NCEL)  for  the  mass  density  of 
each  cell  in  the  HAG  by  appropriately  calling  Subroutine  ECRD  to  get 
successive  columns. 

d.  Changes  the  geographic  location  of  the  event  point  by  a 
small  amount  so  that  if  the  event  point  is  inputted  on  a cell  boundary, 
as  it  would  be  if  it  were  placed,  e.  g.,  at  the  origin  of  the  quadrupole 
coordinate  system,  the  event  point  would  be  moved  off  the  boundary  to 
facilitate  the  code's  assigning  the  event  point  to  the  proper  cell  in  the 
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grid.  This  movement  of  the  event  point  is  made  by  (1)  converting  the 
geographic  coordinates  of  the  event  point  to  quadnapole  coordinates  by 
calling  Subroutine  GEOQUA,  (2)  multiplying  the  quadrupole  coordinates 
by  the  fraction  0.  999997  to  obtain  the  new  quadrupole  coordinates,  and 
(3)  converting  the  new  quadrupole  coordinates  to  geographic  coordinates 
by  calling  Subroutine  QUAGEO. 

e.  Calls  Subroutine  INDEX  to  obtain  the  indices  (lEC,  JEC, 

KEC)  of  the  event  cell  and  the  vertical  fractional-position  parameter 
(FRACB)  for  the  event  point  within  the  event  cell.  Again,  Subroutine 
INDEX  is  different  in  SAIHYD  than  in  NRLHYD.  The  event-point  indices 
are  printed  and  stored  in  EVENTX  Common  as  members  of  arrays. 

f.  Calls  Subroutine  PROPTY  to  obtain,  print,  and  store  in 
array  PP(7)  the  following  event-point  properties:  mass  density,  mass- 
density  scale  height,  heavy-particle  temperature,  and  magnetic-field 
strength,  dip  angle,  declination  angle,  and  dipole  colatitude. 

g.  Moves  the  event  point  slightly  away  from  the  center  or 
boundary  of  a cell  if  either  of  the  quadrupole  coordinates  or  the  geo- 
centric radius  of  the  event  point  differs  from  that  of  the  cell  center  or 

_ 5 

boundary  by  less  than  the  fraction  10  . The  new  coordinate  is  made  to 

_ 5 

differ  from  the  cell- center  or  boundary  coordinate  by  the  fraction  3 10  . 

h.  Computes  the  geocentric  cartesian  coordinates  (XI,  Yl,  Zl) 
of  the  event  point,  given  the  quadrupole  coordinates  and  the  geocentric 
radius  of  the  event  point.  The  Z-axis  is  through  the  quadrupole  north 
pole;  the  X-axis  is  through  the  quadnjpole  west  pole. 

i.  Calls  Subroutine  DEBRIS  which  serves  as  the  principal  con- 
trol routine  for  computing  the  partitioning  and  deposition  of  the  Group-D 
(debris)  energy.  Subroutine  DEBRIS  controls  16  subroutines,  including 
direct  calls  to  Subroutines  HDP.A.RT,  CONSPC,  LEKSPC.  CHXSPC,  and 
BLINE. 
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j.  Sets  up  the  Group-U  (ultraviolet)  properties  by  (1)  calling 
Subroutine  EUXFIT  to  get  the  energy  in  each  of  the  five  Subgroups  L, 

(2)  establishing  pressure  efficiencies  PREFF(L,1)  for  3ubgroup-L 
photons  incident  on  Species  I,  and  (3)  establishing  absorption  coefficients 
XMUl  of  ©2  for  Subgroup  1 of  Group  U.  (These  latter  two  functions 
should  probably  be  done  elsewhere.)  The  calculated  pressure  efficien- 
cies follow  the  somewhat  arbiti'ary  assumption  in  FS-73  that  all  photo- 
electron kinetic  energies  below  10  eV,  and  0.  80  of  the  kinetic  energy  of 
photoelectrons  more  energetic  than  10  eV,  are  thermalized;  the  remain- 
ing energy  is  discarded,  as  having  been  lost  in  radiative  processes. 

k.  Computes  the  deposition  of  the  Group-X  (x  ray),  Group-U, 
and  Group-CHEX  (charge  exchange)  energies  by  making  a triple  DO-loop 
over  all  the  cells  in  the  HAG,  during  which  calls  are  made  to  the  follow'- 
ing  subroutines: 

PLINE  - Computes  line  integrals,  with  aid  of  Subroutines  PINT  and 
HEDGE. 

CHXDEP  - Computes  deposition  of  CHEX  energy. 

HPCHEM  - Computes  prompt  chemistry  due  to  Group-D  energ\' 
deposition. 

PHEAT  - Computes  heating  due  to  Group-D,  -X,  and  -U  energies. 

PCHEM  - Combines  results  of  Subroutines  HPCHEM  with  prompt 
chemistry  due  to  Groups  X and  U. 

l.  Stores  in  array  BUF  those  cell  quantities  updated  in  Sub- 
routine PCHEM  and  stored  in  the  temporary  scratch  area. 

Figure  1 is  a simplified  flow  chart  that  summarizes  the  func- 
tions performed  by  Subroutine  PROMPG.  Figure  2 is  a simplified  flow 
chart  of  the  triple  DO-loop  w'ithin  Subroutine  PROMPG  that  loops  over 
all  cells  in  the  HAG- 
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Compute  event-point  geocentric  cartesian  coordinates 


Call  DEBRIS  for  debris-enorpv  pni-tition  & deposition 


Call  EUXFIT  to  get  Group-U  subgi  '*up  energies 


Set  up  pressure  efficiencies  for  Group-U  subgroups 


Call  PHEAT  to  calculate  specific  energy*  d«»posited  in  event  cell 


r ' 


Triple  DO-loop  to  update  cell  quantities  in  entire  grid 


Move  cell  quantities  from  ten;porary  scratch  area  fo  Time-Slot  2 in  ECS 


, * 

(^RETURN^ 


F’ip.  1.  Simplified  Flow  Chart  for  Subroutine  PROMPG  Which 
Controls  the  Enerpy-Deposition  Routines. 
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Fig.  2.  Simplified  Flow  Chart  for  the  Subroutine  PROMPG 
Triple  DO- Loop  Over  All  Cells  in  the  High-Altitude 
Grid. 
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2.  SOME  AUXILIARY  ROUTINES 


2.  1 SUBROUTINE  BOUNDY 

Subroutine  BOUNDY  establishes  a one-dimensional  array 
GR(NCEL)  of  the  geocentric  radii  to  all  the  cell  top-boundaries  in  the 
HAG.  This  information  is  developed  by  using  the  cell- center  heights 
stored  as  cell  properties  and  obtained  by  calling  Subroutine  ECRD. 

Subroutine  BOUNDY  is  different  in  SAIHYD  and  NRLHYD.  As 
discussed  in  Volume  16,  the  differential-quadrature  formulation 
of  the  Lagrangian  hydrodynamics  (NRLHYD)  does  not  employ  the  concept 
of  a cell  in  computing  the  hydrodynamic  motion,  as  is  done  in  the 
difference-equation  formulation  of  the  Lagrangian  hydrodynamics 
(SAIHYD).  However,  the  deposition  of  energy  (in  SAIROS)  has  been 
formulated  in  terms  of  a cellular  grid  for  use  with  SAIHYD.  To  make 
use  of  this  energy  deposition  procedure  with  NRLHYD,  it  is  necessary 
to  establish  cells  and  cell  boundaries  invNRLHYD  to  correspond  to  the 
same  quantities  in  SAIHYD.  Subroutine  BOUNDY  in  NRLHYD  establishes 
such  cell  boundaries. 

For  N Lagrangian  points  (initially  specified  in  Subroutine 
ATMOSG  to  be  at  the  N-1  SAIHYD  cell  centers  plus  one  more  point  at 
the  SAIHYD  grid  bottom),  the  N-1  interior  boundaries  are  defined  to  be 
midway  between  the  adjacent  Lagrangian  points.  For  the  top  cell,  the 
outer  boundary  is  placed  so  that  the  top  Lagrangian  point  is  at  the  center 
of  the  top  cell.  The  parameter  HABOT,  originally  read  in  by  DRIVER 
in  SETUPl  as  the  bottom  of  the  grid  for  SAIHYD.  is  reset  in  .A.TMOSG  to 
be  the  bottom  boundary  of  the  lowest  NRLHYD  cell.  Since  the  value  of 
HABOT  is  not  allowed  to  change  with  time,  the  bottom  Lagrangian  point 
does  not  remain  at  the  center  of  the  bottom  cell  despite  its  not  moving. 
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SUBROUTINE  INDEX 


2.  2 

For  a given  Point  X in  geographic  coordinates  (OCX,  GLX), 
Subroutine  INDEX  calls  Subroutine  GEOQUA  to  compute  the  quadrupole 
coordinates  (QCX,  QLX)  of  the  point.  Subroutine  INDEX  then  finds  the 
indices  IX,  JX,  and  KX  of  the  Point  X.  If  the  Point  X is  to  the  quadrupole- 
north  of,  quadrupole-west  of,  or  below  the  grid  boundary,  index  IX,  JX, 
or  KX  is  set  equal  to  zero.  If  the  Point  X is  to  the  quadrupole-south  of, 
quadrupole- east  of,  or  above  the  grid  boundaries,  the  index  IX,  JX,  or 
KX  is  set  equal  to  (NCELX+1),  (NCELY+1),  or(NCELZfl),  respectively, 
and  an  appropriate  warning  message  is  printed.  Subroutine  INDEX  also 
computes  the  fractional-height  position  of  the  Point  X within  the  cell  at 
the  time  the  subroutine  is  called. 

If  the  Point  X is  on  a quadrupole  colatitude  or  east  longitude 
boundary.  Subroutine  INDEX  identifies  Point  X as  being  in  the  cell  lor 
which  the  boundary  is  the  south  or  east  boundary  of  the  cell. 

The  version  of  Subroutine  INDEX  used  with  NRLHYD  differs 
from  that  used  with  SAIHYD  mainly  in  that  (1)  the  vertical-direction 
index  (KX)  has  a plus  or  minus  sign  associated  with  it  to  denote  whether 
the  Point  X is  above  or  below  the  Lagrangian  Point  KX,  respectively, 
and  (2)  the  fractional-height  parameter  refers  to  the  fractional  distance 
between  the  two  adjacent  Lagrangian  points  instead  of  the  fractional  dis- 
tance within  the  cell  (with  both  fractional  distances  being  measured 
upward  from  the  bottom). 

2.3  SUBROUTINE  PROPTY 

Subroutine  PROPTY  computes  the  mass  density  [PP(D],  the 
density  scale  height  [PP(2)],  and  the  heavy-particle  temperature  [PP(3)| 
of  a given  Point  X inside  the  HAG.  By  calling  Subroutine  BFIELD. 
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Subroutine  PROPTY  also  obtains  the  local  dip  angle  [PP(4)],  declination  I 

angle  [PP(5)],  strength  [PP(6)j,  and  cosine  of  the  dipole  colatitude 
[PP(7)]  of  the  dipole  magnetic  field  line  passing  through  the  point. 
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3.  DEBRIS  ENERGY  PARTITION  AND 
PARTICLE-SOURCE  SPECTRA 


3,1  INTRODUCTION 

Certain  overall  partitioning  of  the  total  yield  for  three  selected 
devices  is  given  by  Knapp  et  al.  'in  Volume  12 as  i)art  of  the  GET 
ROSCOE  Model  3.  This  partitioning  includes  statements  for  the  frac- 
tional yields  of  x rays  and  hydrodynamics.  The  additit)nal  partitioning 
of  the  hydrodynamics  yield  into  several  types  of  heavy- particle  emis- 
sions and  UV  emission,  needed  for  computing  the  energy  deposition  in 
the  HAG,  is  not  completely  understood.  We  have,  of  necessity,  adopted 
a specific  procedure,  but  it  certainly  deserves  review  and  improvement. 

It  is  convenient  to  classify  the  heavy-particle  emission  from 
the  expanding  weapon  debris-air  mixture  into  loss-cone,  ion-leak,  and 
charge-exchange  (CHEX)  particles  even  though  the  division  among  these 
three  groups  is  not  unique.  For  our  purposes,  we  consider  loss- cone 
particles  to  be  debris  ions  that  escape  the  burst  region  along  the  mag- 
netic field  direction.  Ion- leak  particles  are  also  emitted  along  magnetic 
field  lines  but  they  are  composed  mainly  of  air  ions  heated  by  the  weapon 
debris  expansion.  CHEX  particles  are  mainly  fast  neutral  air  atoms 
born  through  ion-to-neutral  charge- exchange  reactions  within  the  debris- 
air  mixture. 

Partitioning  of  the  weapon  kinetic  yield  into  loss-cone  and  ion- 
leak  energies  follows  the  prescription  given  by  Crevier  and  Kilb  I in  Volume 
10]  as  part  of  the  MRC  ROSCOE  Model  12a.  The  fraction  of  tlie 
remaining  kinetic  yield  emitted  as  CHEX  particles  is  computed  from  the 
RANC-IV  formulation  [GE-70].  Thus,  the  sum  of  the  loss-cone,  ion- 
leak,  and  charge-exchange  yields  is  subtracted  from  the  total 
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hydrodynamic  yield  before  the  UV  yield  is  determined  by  the  prescrip- 
tion of  Fajen  and  Sappenfield  [FS-71a,  FS-73].  We  do  not  discuss  the 
physical  validity  of  our  modeling  of  the  hydro-yield  partitioning;  we  do 
discuss  and  provide  Fortran  energy-deposition  routines.  Hopefully, 
future  modifications  of  the  hydro-yield  partitioning  will  require  only 
minor  modifications  of  these  routines. 

3.2  HYDRO-ENERGY  PARTITIONING 

Partitioning  of  the  weapon  kinetic  yield  is  summarized  in 
Fig.  3;  the  directed  dashed  line  shows  the  computational  order.  The 
energy  indicated  in  each  energy  block  is  successively  subtracted  from 
the  remaining  energy  until  the  residual  UV-plus-thermal  energy  is  ob- 
tained. The  partitioning  into  the  three  heavy-particle  modes  is  per- 
formed by  calls  from  Subroutine  HDPART  to  Subroutines  LOSCON, 
lONLEK,  and  CHXLOS.  The  final  division  into  UV  and  thermal  energies 
is  performed  in  Subroutine  EUXFIT. 


Event 

Kinetic  Yield  | 
I I 


Less-Cone 

Enerpy 


1 J 

1 

Ion- Leak 

CHEX 

|_  1 UV  - Thermal  j 

1 

J 

Energy 

i 

Energy 

1 ! Energy  | 

J 1 ^ 

Fig.  3.  Illustration  of  the  Order  of  Partitioning  of  the  Total 
Kinetic  Yield  into  Loss-Cone,  Ion- Leak,  CHEX, 
and  UV-plus-Thermal  Energies. 
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3.  2.  1 Loss-Cone  Energy 

The  energy  fraction  in  the  loss  cones,  with  neglect  of  collisions, 
is  :Vol.  10,  Eq.  (2-23)] 


Ft,  = 1 - (1  - fj)'  - 12-  (I  - (2  - Cj  -3<^) 


'5e 


(1) 


where 


€ = 


1.40  X 10 


-20 


d a 


K d a 


, :Vol.  10,  Eq.  (2-26):  (2) 


and 


= min(e,  1)  . 


:Vol.  10,  Eq.  (2-22): 


(3) 


Here,  W„  is  the  weapon  kinetic  energy  in  MT,  B is  the  burst-point 

magnetic  field  in  gauss,  v , is  the  debris  velocity  in  cm/sec,  n is  the 

3 Q a 

air-ion  density  in  cm  , Z , is  the  debris-ion  charge  state,  and  Z is 

ct  a 

the  air-ion  charge  state.  However,  with  collisions  included,  the  energ\’ 
in  MT  leaving  the  burst  region  in  the  downward  direction  is  given  by 
"Vol.  10,  p.  14^ 


W 


-i-cd 


--  W F 
2 K ^Ic 


( 

2] 

'n 

max  .0, 

1 - (p  '2  X 10"^^)  : 

\ a ' ' 

1 

L 

(4) 


where  p is  the  burst-point  mass  density,  and  the  energy  leaving  in  the 

cl 

upward  direction  is  given  by  _Vol.  10,  p.  14: 


W,  = 4 F, 
-f-cu  2 K tc 


max 


0,  -il  - (p  2 X 10" 


(5) 


The  total  energy  in  MT  leaving  the  burst  region  through  the  loss  cones 
is  :Vol.  10,  Subroutine  LOSCON  ] 


W 


Ic 


W 


t cd 


W 


-t-cu 


(6) 
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For  more  information  about  this  formulation  the  reader  is  referred  to  . 
Volume  10. 

The  total  fraction  of  the  weapon  kinetic  ener^’  (with  collisions 
included)  leaving  the  burst  retpon  in  the  loss  cones  is 


f,  = W,  /Wj. 
-tc  tc  K 


- 1 2 3 

which,  for  burst-point  densities  less  than  p 2 x 10  s'" cm  , is 

cl 

very  nearly  equal  to  the  fraction  ^ for  neglect  of  collisions.  For  the 
downward  direction  the  fraction  corresponding  to  f^^  is 


f,  . = W,  ,/U'  , . 

Irri  led  K 


For  p < 2 X 10"^^  g/cm^, 

cl 


0.  5 f. 


In  Fig.  4 we  have  plotted  f^^  as  a function  of  the  kinetic  yield  The 

curves  have  been  drawn  through  points  computed  from  the  Subroutine 

LOSCON  given  in  Vol.  10  and  used  in  SAIROS.  The  altitudes  correspond 

to  densities  from  the  ClRA-1965  Model- 5 8-hr  atmosphere  [Cl-65]:  we 

used  B = 0.  5 gauss.  For  the  150-km  curve,  f,  = 0. 

-ted 

3.  2.  2 ton- Leak  Energy 

The  energ>'  in  MT  moving  downward  from  the  burst  region  due 
to  the  ion-leak  mechanism  is  .Vol.  10,  Eq.  (3-1)] 


I I -12 

"'Kd  ' “•‘'^'^1  1’  W?) 


where  Wj,  given  by 


W = W - W 
1 K tc  ’ 


r 


Loss-Cone  Fraction  (f 


k 


is  the  kinetic  yield  left  after  the  loss-cone  energy  is  removed.  No  state- 
ment is  made  in  Vol.  10  about  the  magnitude  of  the  upward-moving  ion- 
leak  energy.  Since  we  must  know  the  amount  lost  by  this  mechanism  in 
all  directions  to  estimate  the  kinetic  yield  left  for  subsequent  yield 
partitioning,  we  arbitrarily  set  the  total  ion-leak  energy  loss  as 

W.,  = 2 W.,  . . (11) 

il  ltd 

Final  determination  of  W.^  must  be  deferred  to  some  future  date. 

The  total  fraction  of  the  weapon  kinetic  energy  leaving  the 
burst  region  by  the  ion-leak  mechanism  is 

f.,  = W.  ,/W^  . (12) 

By  our  simplifying  assumption  in  Eq.  (11),  we  see  that  the  fraction  for 
the  downward  direction  is 

f.,  . - 0.  5 f..  . (13) 

ltd  it 

In  Fig.  5 we  have  plotted  f.^  as  a function  of  the  kinetic  yield  Wj^.  The 
curves  have  been  drawn  through  points  computed  from  the  Subroutine 
lONLEK  given  in  Vol.  10  and  used  in  SAIROS.  We  used  the  same  input 
conditions  as  for  Fig.  4. 

3.2.3  CHEX  Energy 

Time  constraints  did  not  allow  us  to  develop  a new  model  for 
the  energy  leaving  the  burst  region  via  fast  neutrals  produced  through 
charge  exchange.  As  an  interim  procedure  we  have  just  applied  the 
RANG  rv  model  [GE-70]  to  the  remaining  kinetic  yield  so  that  a place 
for  CHEX  particles  might  be  formally  introduced  in  the  code. 
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Kinetic  Yield,  MX 


FiR.  5.  Fraction  of  the  Total  Kinetic  Yield  Escaping  the  Event- Point 
Region  by  the  Ion-Leak  Mechanism,  for  Several  Event 
Altitudes.  Curves  are  drawn  for  B = 0.  5 gauss  and  the 
ClRA-65  mean  atmosphere.  By  our  assumption  the 
downward  fraction  f. , ^ equals  0.  5 f. , . 


The  PIANC  FV''  prescription  for  CHEX-particle  energy  emission 
(in  MT),  adapted  to  our  needs,  is 


W :^-(F  -4F  )W 

^chex  6^CED  ^CEU  ^^CE^'^2 


'^2  = 

where  F„„„,  and  F„„  are  coefficients  associated  with  the 

CrL.1-)  LxliU  Lti 

downward,  upward,  and  horizontal  directions,  respectively.  In  terms 
of  the  burst-point  mass  density  p , their  values  are 
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The  fraction  of  the  weapon  kinetic  energ\-  leaving  the  burst 
region  by  the  charge- exchange  mechanism  is 

f , - W , /W^  . 

chex  chex  K 
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In  Fig.  6 we  have  plotted  a function  of  the  kinetic  yield 

We  question  the  validity  of  the  predicted  CHEX  emission  for  events  be- 
low 150  km.  Some  modification  of  the  model  at  low  altitudes  is  clearly 
required. 


3.2.4  Residual  Energy^  for  UV  and  Thermal  Modes 


After  subtraction  of  the  energies  in  the  loss-cone,  ion-leak, 
and  CHEX  particles,  we  obtain  the  energy  left  for  the  production  of  UV- 
photon  and  thermal  energies.  If  the  residual  energy  is  given  (in  MT) 

by 


= w 
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w.  - w. 


Ic 
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chex 


(17) 


then  the  fraction  of  the  kinetic  yield  available  is 


f 
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W„  W 
K 


K ■ 


(18) 


The  fraction  of  the  weapon  kinetic  energy  available  for  the  UV  and 

thermal  modes  is  plotted  in  Fig.  7 as  a function  of  the  kinetic  vield  W . 

K. 


3.3  PARTICLE  ENERGY  DISTRIBUTION 

Each  source  of  fast  heavy-particles  has  a different  energ>-  dis- 
tribution. We  record  the  distributions  presently  used. 


3.3.1  Loss- Cone  Spectrum 

The  velocity  distribution  of  dowmward-going  loss-cone  parti- 
cles is  assumed  to  be  ^ Vol.  10,  Eq.  (4-4)^ 


dW 

dv 


(MT  sec  cm  ^) 


(19a) 
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Fig.  6.  Fraction  of  Total  Kinetic  Yield  Escaping  the  Event- Point  Regi 
by  the  Charge- Exchange  Mechanism,  for  Several  Event 
Altitudes.  Curves  are  drawn  for  the  same  conditions  as  in 
Phgs.  4 and  5, 


Fi^r.  7.  Fraction  of  the  Total  Kinetic  Yield  Remaining  for  the 
UV'  and  Thermal  Modes  (fR)  After  the  Subtraction  of 
the  Loss-Cone,  Ion- Leak,  and  Charge-Exchange 
Fractions.  Curves  are  drawn  for  the  same  conditions 
as  in  Figs.  4,  5,  and  6. 


Equation  (19a)  is  readily  integrated  to  give  an  integral  velocity  distribu- 
tion which  may  be  transformed  to  an  integral  energy  distribution. 
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j = max,0,  v^  - 4 X 10  , (cm/sec)  iVol.  10,  Eq.  (4-3)  (20d) 
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(20e) 


r 


\'2  = f 2 X (cm/sec)  'Vol.  10,  Eq.  (4-3)] 

8 0 27 

- 2.7X10  (cm/sec),  .Vol.  10,  Eq.  (4-1)]  (20f) 

with  the  debris-ion  mass  in  grams  and  the  burst  kinetic  energy 
in  MT. 


3.3.2  Ion- Leak  Spectrum 

The  velocity  distribution  of  downward-going  ion- leak  particles 
is  assumed  to  be  ]Vol.  10,  Eq.  (4-11)] 
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dv 
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(21) 


Equation  (21)  is  readily  integrated  to  give  an  integral  velocity  distribu- 
tion which  may  be  transformed  to  an  integral  energy  distribution. 
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(23c) 

with  m 

the  ion- leak  particle  mass,  v.^  given  by  Eq.  (20e),  and 

v„  = l.Sxlo"^  (cm/sec)  . ]V^ol.  10,  Eq.  (4-11)' 

(24) 
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3.3.3  CHEX  Spectrum 


The  charge- exchange  particle  spectrum  tends  to  vary  widelv 
with  burst  altitude,  direction  of  emission,  and  the  weapon  expansion 
velocity.  To  satisfy  our  immediate  needs  we  have  assumed  an  expo- 
nential energy  distribution.  Examination  of  earlier  work  [HL-71, 
HM-72b]  tends  to  lend  support  for  this  functional  form.  Explicitly,  the 
integral  energy  spectrum  is 
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corresponding  to  a mass- 15  particle,  and  \'2 


given  by  Eq.  (20e). 
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4.  UV  SOURCE  CHARACTERISTICS 


We  have  adopted  the  HARC-code  prescription,  piven  in  our 
Subroutine  EUXFIT,  for  the  UV  spectrum  as  a function  of  the  event-point 
density  and  kinetic  energy.  This  prescription,  based  on  calculations  by 
Fajen  and  Sappenfield  [FS-71a,  Table  IV:  FS-73  |,  gives  the  fraction  f. 

(i  = 1,  5)  of  the  kinetic  energy  converted  to  UV  in  each  of  five  spectral 
energy  groups  (with  energy-bin  edges  at  0,  13.  6,  15.6,  17,  and  30  eV), 
according  to  the  expression 


f. 
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b.x 
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e.x  + f.)y 


h.x 

1 


(26) 


where 


X = In  K 


y = In  p 

K (Kl’)  = available  kinetic  energy 

3 

p(g/cm  ) = event-point  density  . 


Because  of  the  approximate  nature  of  the  equations,  the  total  fraction 
(FRU)  is  given  [FS-71a]  by  Eq.  (26)  with  another  set  of  coefficients 
(i  = 7 in  the  code).  However,  the  code  limits  FRU  to  a maximum  value 
of  0.98.  If  either  f.  or  FRU  is  negative  for  certain  conditions,  it  is 
set  equal  to  zero.  The  intended  [FS-71al  range  of  validity  of  Eq.  (26) 
is  given  by  the  expressions 

8 X 10“^^  s.  p 6 X 10“^'^ 

5 ^ K s lo'^  . 
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In  the  HARC  Subroutine  DSFIT,  which  we  renamed  EUXFIT 
after  making  minor  changes,  D.  H.  Sowle  introduced  a sixth  energy 
group  due  to  effective,  low-energy  x rays.  The  fraction  fg  is  given  by 

,,  FRX  X EX 
^6  ^6  ^ EUX 


where 


f ' 


0.01 


0.  01  + 


FRU  X EH 
EUX 


f. 


1=1 


EH  = K,  EX  is  the  x-ray  energy,  and  EUX  is  the  sum  of  the  Group- U 
and  an  effective  low-energy  portion  of  Group  X,  given  by 

EUX  = FRU  X EH  + FRX  X EX  . 

Here,  FRX  is  given  by 

FRX  = Max[0.  08(1  - TBAREX),  10"®  | . 

TBAREX  = fraction  of  x-ray  energv’  escaping  a 10-km  radius,  obtained 
by  calling  Subroutine  XTCOEF. 

The  fractions  f.  in  Eq.  (26)  are  normalized  bv  the  expression 

,,  _ h FRU  X EH  , = 

1 5 EUX  ’ 1 1,3  . 

0.01  , ^ f. 

i = l 


One  can  verify  that 

1 - 1 

The  Fortran  variables  for  f ' (i  = 1.  5)  and  fl  are  FR(I),  I = 1.6.  F'R(6) 

i 6 

is  not  used  in  the  HAG  module. 


5.  X-RAY  SOURCE  CHARACTERISTICS 
AND  ENERGY  DEPOSITION 

To  facilitate  development  of  the  HAG  module,  we  used  the 
HARC  Subroutine  XTCOEF  instead  of  the  corresponding  GET  routines 
for  the  x-ray  source  and  deposition  characteristics.  By  improving  the 
HARC  procedure  for  x-ray  energy  deposition,  we  also  eliminated  use 
of  the  HARC  Subroutine  XDCOEF,  used  for  x-ray  energ\'  deposition. 
Our  calculation  of  the  x-ray  energy  deposition  is  done  in  Subroutine 
PHEAT. 

The  HARC  Subroutine  XTCOEF(EM,  TBARX)  provides  the  ef- 
fective energy  transport  coefficient  for  six  spectral  tvpes,  each  char- 
acterized by  the  value  of  an  index  (1  ^ IDSX  • G).  The  arg-ument  EM 
2 

(g/cm  ) is  the  intervening  air  mass  and  TBARX  is  the  fraction  of  the 
energy  reaching  the  depth  EM.  To  prevent  Subroutine  XTCOEF  from 
returning  a physically  impossible  value  in  excess  of  1.0  (due  to  using 
the  same  formula  for  extrapolation  as  interpolation),  we  revised  the 
routine  by  using  the  formula 

F\T 

TBARX  - 1 - [1  - COEFd,  IDSX)] 

10"^ 

for  values  of  EM  less  than  10"'^  g/cm^,  where  COEF(l,  IDSX)  is  the 
coefficient  for  the  spectral  index  IDSX. 

To  deposit  x-ray  energy,  we  use  the  formula 


DAIXX(erg/g) 


EX  TBXIN  - TBXOUT 


(p  OlfR*"  , - IC  ) 
out  in 


where  EX  is  the  total  x-ray  energ\’,  TBXIN  and  TBXOUT  are  obtained 
by  calls  to  Subroutine  XTCOEF(EMIN,  TBXIN)  and  XTCOEF(EMOUT, 
TBXOUT),  and  the  input  arguments  EMIN  and  EMOUT  are  given  bv 
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EMIN(g/cm^)  = f p dR 

■'o 


EMOUT(g/cm^)  = f P dR  , 

•'o 

3 

with  p(g/cm  ) being  the  total  mass  density.  R.  and  R ^ are  the  dis- 

in  out 

tances  from  the  event  point  to  the  points  of  entry  and  exit  of  a cell, 
measured  along  the  ray  path  from  the  event  point  through  the  cell  center. 

For  permanent  use  one  needs  to  replace  the  HARC  Subroutine 

XTCOEF  with  the  GET  " Vol.  12^  equivalent  Subroutines  WOXl(l)  and 

WOXCd,  PMASS,  FCONT)  and  Function  EXP1KTP(X1,  Yl,  X2,  Y2,  X). 

Our  DRIVER  should  call  WOXl(I)  to  read  (a)  the  fractional  x-ray  energy 

in  each  of  18  energy  groups  (SPECX),  and  (b)  the  yield  fraction  in  x rays 

(FEDX),  and  fill  WOX  Common.  One  needs  to  change  the  current  calls 

to  Subroutine  XTCOEF  in  Subroutines  EUXFIT  and  PHEAT  by  calls  to 

Subroutine  WOXC,  and  take  account  of  the  fact  that  FCONT  returned  by 

Subroutine  WOXC  is  the  fraction  of  x-ray  energy  contained  within 
2 

PMASS  (g/cm  ) rather  than  being  the  fraction  transmitted  beyond  EM 
2 

(g/cm  ) as  in  Subroutine  XTCOEF.  One  also  needs  to  get  the  total  en- 
ergy WB  (MT)  into  DEVICE  Common,  as  well  as  convert  from  MT  to 
ergs. 
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6.  DEPOSITION  OF  DEBRIS  PARTICLES 
6.  1 INTRODUCTION 

Before  we  proceed  with  the  detailed  physical  and  geometric 
considerations  necessary  for  deposition,  we  will  outline  the  overall  pro- 
cedure. The  total  weapon  kinetic  yield  is  divided  into  loss-cone,  ion- 
leak,  CHEX,  UV,  and  thermal  energies.  Each  of  these  energy 
components,  except  the  last,  is  deposited  separately  in  the  HAG.  First, 
the  loss-cone  particles  are  deposited  in  a field-aligned  tube  centered 
about  the  burst-point  field  line.  Next,  the  ion-leak  particles  are  also 
deposited  within  a field-aligned  tube  l)ut  of  size  different  from  the  loss- 
cone  tube.  Then  we  deposit  the  CHEX  particles  as  neutral  particles 
without  any  consideration  of  subsequent  ionization  and  magnetic  field 
focusing.  Deposition  of  UV  energy  is  discussed  in  Section  9. 

Changes  of  chemical  species  produced  in  each  cell  by  the  sum 
of  the  heavy-particle  energ>'  depositions  modifv  the  species  composition 
that  would  normally  be  seen  by  the  UV  photon  flux.  Although  the  final 
species  complement  may  differ  in  a heavily-dosed  cell  if  the  UV  were 
deposited  first,  the  differences  are  not  expected  to  be  large.  In  lightly- 
dosed  cells,  the  order  of  irradiation  is  probably  of  minor  consequence. 
At  any  rate,  we  found  it  convenient  to  deposit  the  heavy-particle  energ\- 
first  in  our  implementation. 

6.2  HEAVY- PARTICLE  ENERGY  LOSS 

According  to  heavy-particle  stopping  theory  [LS-63].  the  dif- 
ferential energy  loss  rate  of  a particle  traversing  matter  is 

dE 

-j—  = - n S(E)  (eV  cm)  i27) 
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where  E is  the  particle  energy,  x the  distance,  and  n the  number  den- 
sity of  target  atoms.  The  stopping  cross-section  S(E)  has  an  elastic 
contribution  ol  and  an  inelastic  portion  viz, 


S(E)  = a + 


(eV  cm' 


(28) 


Values  for  oxygen  and  aluminum  atoms  slowing  down  in  nitrogen  are 
found  by  direct  substitution  into  Eqs.  (2.  5)  and  (2.  7')  in  LS-63.  If  the 
particle  energy  E is  measured  in  eV,  the  values  of  a and  ^ are,  for 
oxygen. 


= 2.  98  X 10 


[eV  cm^  (target  atom)"^ 


= 1.  13  X 10 


-16 


t 2 “1 

[eV  cm  (target  atom) 


and,  for  aluminum. 


“ai  ^ '^0 


-14 


2 - 1 

[eV  cm  (target  atom) 


= 1.  92  X 10' 


[eV*"  cm^  (target  atom)  ^ 


Note  that  the  value  of  is  much  larger  than  that  propagated  into  the 
work  of  HM-72b  from  Appendix  4A  of  HL-71. 

In  Eq.  (27)  it  is  often  convenient  to  introduce  the  areal  mass 
2 

density  (g/cm  ) 

d4  = mn  dx  (29) 

traversed  by  a particle  in  lieu  of  the  rectilinear  distance;  i.  e.  . 

(30) 


m ^ = - S(E)  . 


Here,  m is  the  mass  of  a target  atom.  From  Eq.  (30)  we  can  compute 
by  straightforward  integration  the  areal  mass  densitv  required  to  slow 
a particle  from  energv-  E^  to  E2. 
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Ola) 


After  inserting  Eq.  (28)  in  Eq.  (31a)  and  integi’ating,  we  obtain  the 
result 


(31b) 


which  is  required  for  the  calculation  of  energy  deposition  in  a cell. 

In  our  application  of  Eq.  (31b)  we  will  normally  know  ^ along 
the  particle  path  in  a cell  and  the  cell-entrance  energy  E^^  of  the  parti- 
cle; we  will  then  calculate  the  particle  cell-exit  energ>^  by  use  of 
the  Function  Routine  FZET  which,  given  the  range,  iteratively  solves 
for  the  energy  in  the  range- energy  equation.  The  difference  E^  - E^ 
is  then  the  total  energy  deposited  by  the  particle  in  the  cell,  of  which  a 

fraction  f is  due  to  elastic  collisions  and  a fraction  f „ is  due  to  in- 

a p 

elastic  collisions,  where 


f 

a 


5E^/(E^-E2) 


(32a) 


6E  = a 4(E.  E„)/m 

a 1 2 


(32b) 


f 


)3 


1 - f . 


a 


(32c) 


If  the  error  in  using  the  first-order  Taylor  e.xpansion  of  Eq.  (29) 
is  less  than  the  corresponding  tolerance  accepted  in  using  the  Function 
Routine  FZET,  as  is  the  case  if  the  traversed  areal  density  is  verv 
small,  we  compute  the  cell-exit  energy  by  using  the  (second-order) 
Taylor  expansion  of  Eq.  (30). 
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6.  3 ENERGY  DEPOSITION  OVER  A DISTRIBUTION  FUNCTION 


Although  Eq.  (28)  gives  the  stopping  cross-section  for  a mono- 
energetic  particle  of  energy  E,  the  particle  emissions  span  a distribution 
of  energies.  The  sum  of  the  elastic  (e^)  and  inelastic  (e.)  energ\'  con- 
tributions per  target  atom  (eV  atom"^)  is 


e 

t 


CO 


dE  \ 
dN  o| 

dE^  dE  L ^ 
o /E  =E  -E 


(33a) 


(33b) 


with  A the  cross-sectional  area  of  the  particle  beam  at  the  deposition 

location,  E^  the  particle  energy  at  the  source  prior  to  penetration  of  any 

matter,  E the  particle  energy  at  the  deposition  location  after  traversal 
2 

of  ^(g/cm  ) of  material,  and  E^  the  energy  lost  by  the  particle  to  reach 
the  deposition  location.  The  energy  transformation  function  is  just 

dE  S(E  ) 
dE  “ 'sIeT 

and  the  number-vs-energy  distribution  is  related  to  the  integral  number 
spectrum  by 

00 

= -i  w-  ■ 

E o 


6.4  PRACTICAL  IMPLEMENTATION 

In  the  practical  situation  where  deposition  must  be  done  in 
hundreds  of  cells,  it  is  desirable  to  have  analytical  expressions  for  the 
two  integrals  in  Eq.  (33a).  Number-distribution  functions  for  the  particle 
sources  outlined  in  Section  3.3  do  not  lead  to  expressions  analytically 
integrable. 
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To  circumvent  this  difficulty  we  have  represented  each  con- 
tinuous energy  distribution  by  a number-vs- energy  histogram  of  equal 
energy  intervals.  Each  bar  of  the  histogram  represents  a fixed  number 
of  particles  at  a single  energy  midway  between  the  edge  values.  Each 
monoenergetic  beam  is  then  deposited  as  a separate  source.  An 
attractive  feature  of  this  approach  is  its  adaptability  to  any  distribution 
that  may  arise  in  the  future. 


6.4.1  Histogram  Representation 

If  the  number-vs- energy  distribution  of  the  source  function  is 
represented  by  p monoenergetic  beams, 

f = X;  N,  S(E  - E^)  , (36) 

k=l 

then  the  elastic  and  inelastic  energies  (eV  atom"^)  deposited  at  a depth 
2 

^(g/cm  ) from  the  source  location  are,  from  Eq.  (33b), 


and 


k=l 


P 


k=l 


00  P 

i = IX  ® ® E N, S(E  - 

k=l 


lE  - : 

k-1 


(37a) 


(37b) 


(38a) 


(38b) 


L 
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The  number  of  particles  at  energy 
is  found  from 


for  each  particle  source  function 


N 


k 


i yL 


(39) 


where  W assumes,  in  turn,  each  of  the  values  (expressed  in  eV)  of 

W,  ,,  V/. , and  W . . Beam  enere;ies  E,  are  determined  from  the 

-ted  i-td  chex  k 

p+1  histogram-bar  edges  by 


E 


k 


i = 1,  . . . , p 
k = i . 


(40) 


6.4.2  Loss- Cone  Histogram 

Histogram-bar  energy  edges  for  the  downward- going  loss-cone 
yield  is  found  by  equating  the  right-hand  side  of  Eq.  (19b)  to  the  fraction 
of  the  energy  below  Explicitly  we  have 


e5/2  _ ^5/2 
1+1  min 
e5/2  _ g5/2 
max  min 


1 

P 


i = 1,  . . . , p 


(41) 


or 


'i*l 


with 


.i( 

„5/2 

E 

2/5 

min  p \ 

max 

min/ 

E . , 

E = E 

min 

P+1 

max 

i = 1,  . . . , p (42) 


6.4.3  Ion- Leak  Histogram 

For  the  ion- leak  source  we  use  Eq.  (22)  to  find  the  energ\'  edges. 
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with 


E.^  1 - E^  . 

1 ' 1 min 

E^  -E^ 
max  min 


1 

P 


i = 1,  . . . , p 


(43) 


or 


E.  = 

r 

*■ 

E‘  . 

— 1 

1 i 

;e" 

E*.  )] 

1 

min 

P ' 

^ max 

min^ 

i = 1,  . . . , P (44) 


E,  = E . , 

1 nun 


E , = E 
p-1  max 


6.4.4  CHEX  Histogram 

In  a similar  manner  we  use  Eq.  (25a)  to  find  the  energy 
histogram-bar  edges  for  the  CHEX-particle  source: 


-aE  -aE,  , 

„ min  1-1 

J ^_e 

-aE  -aE 

^ min  max 

e - e 


i_ 

P 


i = 1, 


(45) 


or 


with 


E.  . = - - In 

1+1  a 


-aE  . I -aE  -aE 

min  1 min  max 

e — (e  - e 

P \ 


i = 1,  . . . , p 


(46) 


E,  = E . 


min 


E 1 = E 
p+l  max 


6.  5 PARTICLE  SPATIAL  DISTRIBUTIONS 

In  addition  to  an  energ\'  distribution,  each  of  the  three  heavv- 
particle  sources  has  its  own  spatial  distribution.  Since  we  did  not  have 
time  to  introduce  the  re-direction  effects  of  the  geomagnetic  field  in  the 
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deposition  of  CHEX  particles,  we  assumed  the  particles  to  be  emitted 
with  spherical  symmetry.  However,  the  loss-cone  and  ion-leak  parti- 
cles are  guided  by  the  geomagnetic  field  in  the  deposition  processes 
within  tubes  about  the  burst-point  field  line.  In  this  section  we  discuss 
the  deposition  and  the  particle  spatial  distribution  within  these  tubes; 
we  ignore  spiralling  effects. 

Before  presenting  the  formulas,  we  describe  the  overall  pro- 
cedure. The  differential  foi'm  for  the  radial  distribution  of  the  downward- 
moving  loss-cone  and  ion-leak  energies  has  been  specified  at  100-km 
altitude  by  Crevier  and  Kilb  Vol.  10".  We  consider  the  corresponding 
integral  distribution  function  and  form  q concentric  rings  in  the  tangent 
plane  so  that  each  ring  contains  1 q of  the  total  energ\’  for  either  the 
loss-cone  or  ion-leak  particles.  We  compute  the  location  (geographic 
colatitude  and  east  longitude)  of  the  center  of  the  distribution  function, 
i.e.  , the  origin  of  the  tangent-plane  coordinate  system  at  100-km  alti- 
tude, by  tracing  the  event-point  field  line.  The  source  points  (in  the 
vicinity  of  the  event  point)  corresponding  to  representative  points  in  the 
concentric  rings  are  determined  by  tracing  the  field  lines  passing  through 
the  representative  points  to  their  intersections  with  the  equipotential 
surface  through  the  event  point. 


6.  5.  1 Representation  of  the  Spatial  Distribution  at  100-km  Altitude 

The  spatial  distribution  suggested  by  Vol.  10  for  the  downward- 
going particle- energy  in  the  loss-cone  and  ion- leak  tubes  is  of  the  form 


= :4 


ttR 


o L 


1 t (r  R^)2 


dr  do  , 


(47) 


such  that  the  integration  over  the  cross-sectional  area  of  the  tube. 
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(48) 


W 


7TR‘'  *'0 

o 


/ ‘^  y 

-'A  •'n 


dr 


(r/R  )2 
o 


= W 


is  just  the  downward- going  energy.  In  these  expressions  r and  <p  are 
the  coordinates  in  the  tube's  cr'^ss-sentional  area,  acquires  the 
value  or  and  R^  acquires  the  characteristic  tube  radius 

R^^  or  R.^  for  the  loss-cone  and  ion-leak  particles,  respectively.  The 
variable  r is  ultimately  defined  in  Vol.  10  as  the  horizontal  distance,  at 
100-km  altitude,  from  the  burst-point  field  line  to  the  field  line  through 
a representative  point. 

The  integral  distribution  function  for  the  particle  energ>'  passing 
within  a distance  r of  the  event-point  field  line  is  found  from  Eq.  (48)  by 
integrating  from  0 to  radius  r,  i.  e.  , 


W^(0,r) 


W 


= 1 


1 


d 


r2,R2 

0 


(49) 


To  avoid  the  concept  of  a tube  with  infinite  radius,  we  limit  the  radius 
of  the  tube  to  3R^.  According  to  Eq.  (49),  W^(0,  3R^)/W^  = 0.  90,  i.e. , 
0,  90  of  the  energy  lies  within  the  radius  3R^.  However,  we  will  re- 
normalize the  distribution  function  so  that  the  entire  energ>’  is 
within  the  radius  3R^.  Our  new  integral  distribution  function  for  the 
radial  distribution  of  the  particle  energy  is 


f = 


(50al 


where 


(50b) 
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If  we  solve  Eq.  (50a)  for  p, 


0.  9f  _ 
1 - 0.  9f 


(51a) 


or 


0.  9f 
1 - 0.  9f 


(51b) 


we  then  can  compute  the  radius  for  a given  value  of  f.  For  the  present 
we  are  using  four  annular  rings,  with  outer  radii  corresponding  to 
values  of  f equal  to  1/4,  1/2,  3/4,  and  1,  and  characteristic  circles 
within  these  annular  rings  with  radii  corresponding  to  values  of  f equal 
to  1/8,  3/8,  5/8,  and  7/8.  These  radii  are  tabulated  in  Table  1.  On 
each  characteristic  circle  we  place  four  representative  points  at  90-deg 
intervals  except  that  we  also  place  a representative  point  at  the  center 
and  on  the  fourth  characteristic  circle  we  double  the  number  of  points. 
To  the  set  of  s equally- spaced  representative  points  on  alternate  char- 
acteristic circles,  we  apply  a rotation  about  the  event-point  field  line 
equal  to  n/s.  The  weight  assigned  to  each  of  the  s representative 
points  within  one  of  the  q annular  rings  is  l/(qs).  The  annular  rings, 
characteristic  circles,  and  representative  points  currently  used  are 
shown  in  Fig.  8. 


6.  5.  2 Particle  Allocations 

The  spat’al  allocation  of  the  source  particles  (cf.  Eq.  (39)) 
at  energy  Ej^  is  simply  Nj./ (qs)  for  every  field  line  traced. 
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Table  1.  Radii  of  Annular  Rings  and  Characteristic  Circles  Used 
to  Approximate  the  Radial  Distribution  F'unction  of 
Downward-Moving  Loss-Cone  and  Ion- Leak  Energies. 


Integral  Fraction 
of  Energy,  f 

Outer  Radius  of 
Annular  Ring^ 

Radius  of 

Characteristic  Circle^ 

0 

0 

0. 125 

0.  356 

0.250 

0.  539 

0.375 

0.  713 

0.  500 

0.  904 

0.  625 

1.  134 

0.  750 

1.441 

0.  875 

1. 925 

1.000 

3.  000 

In  units  of  the  characteristic  tube  radius  R^. 


6.6  STARTING  COORDINATES  OF  THE  FIELD  LINES 

It  remains  to  establish  the  geocentric  coordinates  (longitude, 
colatitude,  and  geocentric- radius)  of  each  of  the  particle  source  points 
in  the  vicinity  of  the  event  point.  As  mentioned  above,  these  source 
points,  corresponding  to  the  representative  points  in  the  tangent  plane 
at  100-km  altitude,  are  determined  by  tracing  the  field  lines  passing 
through  the  representative  points  to  their  intersection  with  the  equipo- 
tential  surface  through  the  event  point.  We  now  derive  the  equations  to 
determine  these  intersections. 

The  equation  for  a dipole  field  line  is 

2 

r = r sin  \ (52) 

o 


42 


V,  North 


P’lft.  8.  Four  Annular  Kin^s,  Four  Characteristic  Ci  rck’s,  and  21 
Representative  Points  Used  to  Approximate  the  lUidial 
Distribution  of  Dounward-Movinii  Loss-Cone  and  Ion- Leak 
Fner^ies.  The  plane  of  the  paper  represents  a tanuent  plane 
at  100-km  altitude  with  orittin  at  the  event-pomt  fiidd  line. 
Fast  on  the  diaitrani  corresponds  to  the  tanuent-plane 
X axis. 
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where  is  the  equatorial  radius  and  \ is  the  dipole  colatitude.  (The 


r in  Eq.  (52)  should  not  be  confused  with  the  r in  Section  6.  5.  1.)  The 
potential  function  for  a magnetic  dipole  is 


V = ' M ( r cos  \ 


where  |m!  is  the  magnitude  of  the  magnetic  dipole  moment.  For  a 


constant  value  of  V, 


- 2 - 2 
Ml  r cos  \ = 1m1  r^  cos  , 


where  the  subscript  1 denotes  the  event  point,  so  that  the  equation  for 
the  equipotential  curve  is 


cos  \ 


1 cos  X, i 

L ij 


The  equation  for  the  field  line  through  a source  point  is  given 
by  Eq.  (52)  which  we  rewrite  as 


cos  X = 1 - I’/i’ 

o 


By  eliminating  cos  \ between  Eqs.  (55)  and  (56),  we  get 
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B = 1 + h /R  . (59c) 
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We  use  Newton's  method  of  iteration  to  solve  the  quartic,  Eq.  (58),  for 
the  (reduced)  geocentric  radius  at  which  the  source-point  field  line 
intersects  the  event-point  equipotential  sui'face.  Given  the  geographic 
coordinates  of  a representative  point  in  the  tangent  plane  at  100-km 
altitude,  Subroutine  CONJUG  will  trace  the  field  line  through  the  point 
to  an  altitude  h^  and  return  the  desired  geocentric  coordinates  (latitude 
and  longitude)  of  the  intersection  of  the  source-point  field  line  and  the 
event-point  equipotential  surface. 

6.7  FIELD- LINE  tracing 

We  need  to  know  the  heavy-particle  energy  deposited  in  each 
cell,  which  can  be  found  if  we  know  the  areal  mass  density  along  the 
path  within  a cell  (the  product  of  the  field- line  path  length  within  the  cell 
and  the  mass  density  in  the  cell).  Thus,  we  must  trace  a field  line 
from  the  source  point  through  each  cell  until  the  boundary  of  the  HAG 
is  reached. 

To  trace  the  field  lines  we  start  at  a given  source  point  at  alti- 
tude hg.  Given  the  geographic  coordinates  of  the  source  point  and  a 
new  altitude  h_  and  altitude  step  Ah,  where 

V 

h = h - Ah  , (60) 

vs 

Ah  = At.  sin  I , (61  i 

and  I is  the  local  dip  angle.  Subroutine  CONJUG  will  return  the  geo- 
graphic coordinates  of  the  new  altitude  point  and  Subroutine  INDE.X  will 
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return  the  indices  of  the  cell  containing  the  new  point.  If  the  new  point 
is  in  the  same  cell  as  the  source  point,  a second  new  altitude  point, 
lower  by  an  amount  <^h,  is  selected.  Again,  Subroutine  CONJUG,  given 
the  geographic  coordinates  of  the  source  point  and  the  altitude  fh^  - 2Ah) 
of  the  second  point,  will  return  the  geographic  coordinates  of  the  second 
point  on  the  field  line  and  Subroutine  INDEX  returns  the  indices  of  the 
cell  containing  the  second  point.  This  procedure  is  continued  until  the 
cell  indices  indicate  the  field  line  has  exited  from  the  cell.  An  iterative 
backup-and-boundary-crossing  procedure  is  then  followed  in  which  A-t, 
initially  set  to  8 km,  is  repeatedly  halved  and  the  boundary  crossing  is 
finally  accepted  when  At  = l km.  The  field- line  path  length  from  the 
starting  point  to  the  (approximate)  exit  point  is  also  returned  by  Sub- 
routine CONJUG.  This  exit  point  is  then  used  as  a new  starting  point 
with  At,  = 8 km  and  the  field  line  is  further  traced  until  it  exits  the  cell 
from  which  it  last  started.  This  procedure  is  continued  until  the  field 
line  leaves  the  HAG  or  until  the  most  energetic  particle  being  followed 
has  stopped.  This  procedure  is  then  repeated  for  everv  source  point 
corresponding  to  the  set  of  representative  points  in  the  tangent  plane  at 
100- km  altitude. 


6.  8 LOSS-CONE  AND  ION- LEAK  PARTICLE  DEPOSITION 


Energv  deposition  for  loss-cone  and  ion-leak  iiarticles,  while 
based  on  Eqs.  (371i)  and  (38b),  must  be  done  different Iv  to  conserve 
energy  explicitly.  In  particular,  the  specific  energies  (erg  g)  deposited 
in  a cell  bv  loss-cone  and  ion-leak  heavv-particle  elastic  collisions  and 
inelastic  collisions,  are,  respectivelv, 
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and 
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= aF  y NJE.  - E (1  - f 
1 M » k in  out  k c 

k=l 


where  E.  and  E ^ are  the  entrance  and  exit  energies  (eV)  of  the  par- 
in  out 

tide  in  the  kth  energy  group,  M is  the  cell  mass  in  grams, 

= 1.  6 X 10  erg/eV,  and  f^j^  is  given  by  Eq.  (32a)  with  an  addi- 
tional subscript  k. 

6.9  CHEX- PARTICLE  DEPOSITION 

As  mentioned  in  Section  6.  5,  the  CHEX-particles  are  considered 
emitted  with  spherical  symmetry,  or  the  energy  emitted  per  steradian 
is  given  by  ^'chex'  Particle  paths  through  each  cell  are  straight  and 
have  a direction  determined  by  the  line  that  connects  the  event  point  and 
the  center  of  the  cell.  Path  lengths  within  each  cell,  from  entrance  to 
exit,  are  found  by  the  Subroutine  PLINE  and  the  intercepted  steradianal 
mass  is  given  by 

M„  = ip(R^  . - ) • (64) 

D 3 out  m 

Here,  p is  the  cell  mass  densitv;  R.  and  R ^ are  slant  distances  from 
’ in  out 

the  event  point  to  the  cell  entrance  and  exit,  respectively. 

The  specific  energies  deposited  in  a cell  by  CHEX  heavy- 
particle  elastic  and  inelastic  collisions,  are,  respectively, 

e P N 

e'  = ^ V ^ (E.  - E ),  f , (65) 

e M,^  ^ 4tt  in  out  k ok 
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r 

I where,  again,  E.^  and  are  the  entrance  and  exit  energies  (eV)  of 

I the  particle  in  the  kth  energy  group. 


i 

» 
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7.  PARTITIONING  OF  THE  DEPOSITED  DEBRIS  ENERGY 
AND  CONSEQUENT  SPECIES  (SUBROUTINE  HPCHEM) 

Subroutine  HPCHEM  partitions  the  total  ener^  lost  in  each  cell 
by  loss-cone,  ion- leak,  and  CHEX  particles  undergoing  inelastic  and 
elastic  collisions.  The  energy  is  partitioned  into  heavy-particle  thermal, 
electron  thermal,  dissociation,  ionization,  excitation,  and  radiation.  A 
model  of  heavy- particle  prompt  chemistry  is  also  used  to  compute  the 
changes  in  number  densities  of  N2,  02,  N(4S),  N{2D),  O,  N2+,  02+,  N*, 
0+,  and  electrons. 

One  of  the  simplifications  made  in  depositing  the  heavy-particle 
energies  is  that  the  range- energy  relation  used  in  Subroutines  DEPO 
and  CHXDEP  is  based  on  the  total  mass  density  of  a cell  whereas  in 
Subroutine  HPCHEM  only  the  species  N2,  02,  and  O are  considered  as 
targets.  The  modeling  of  the  associated  heavy- particle  prompt  chem- 
istry is  based  on  a speculative  treatment  [HM-72b.  Section  6:  HL-73b, 
Section  5],  modified  by  the  need  to  interface  with  the  late-time  grid 
chemistry.  (The  rationale  for  the  heavy-particle  prompt-chemistry 

modeling,  to  the  extent  it  exists,  wall  not  be  repeated  here.)  In  this 

2 

late-time  chemistry,  for  atomic  nitrogen  only  N(  D)  is  carried  as  an 
excited  term  and  for  atomic  oxygen  only  0(^D)  is  accounted  for  (im- 
plicitly)  as  an  excited  term.  Thus,  N(  P)  and  0(  S)  are  purposely  ex- 
cluded as  possible  reaction  products  of  the  heavy- particle  prompt 
chemistry.  Conversion  of  N^  and  into  NO^  (the  only  molecular  ion 
in  the  grid  chemistry)  is  performed  at  the  beginning  of  Subroutine 
PCHEM. 

In  Table  2 we  summarize  the  reactions  from  inelastic  and 
elastic  collisions  that  are  modeled  in  the  prompt  cliemistry  from  heavy- 
particle  energy  deposition. 
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Table  2.  Heavy- Particle  Prompt-Chemistry  Reactions 
from  Inelastic  and  Elastic  Collisions. 


Number  of 

Branching 

Reactions,^ 

Reaction 

Ratio 

cm"  3 

1 

i 

N^.e 

1/1.35  j 

|n(^s)  + 

+ e 

0.35/1.35  ) 

e 

n ^ 

y ' 

H 

OT  , 

i ^ 

1 

O^.e 

1 1.35  ) 

2[0  ] 

w 

) 3 

1 

|°2" 

[OCP)  + 

+ e 

0.35/1.35  ) 

o - 

\ 

O + e 

i^tel 

n ' ‘ 

1 

4 

( N(^S) 

4 14  j 

2[N2] 

NfS)  - 

) 9 

( 

n 

n c 

1 

( N(^D) 

10/14  ) 

1 

o 

H 

(0(^P) 

9/14  j 

2[0  ] 

m ( 

|°2- 

orp)  + 

I 0(^D) 

( 

n 

n c 

1-^ 

u 

5/14  ) 

1°  ■ 

j 0(^P) 

9 14  j 

1 0(^D) 

5/14  1 

n "c 

^ n = 2|[N2]  ' l02](  + [O]  = number  density  of  nuclei, 
n^,  = number  density  of  particles  in  cascade  of  elastic-cnllisi^ni  enercy. 
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In  Subroutine  HPCHEM,  before  dealing  explicitly  with  the  par- 
titioning of  the  species  and  energies  from  the  inelastic  and  elastic  colli- 
sions, we  must  make  a number  of  preliminary  steps  and  establish 
several  parameters  used  in  the  modeling. 

1.  From  the  extended  core  storage  for  the  cell  of  interest, 

3 -3 

obtain  the  total  mass  density  P(g/cm  ) and  the  number  densities  (cm  ) 

of  N2,  O^,  and  O. 

2.  For  both  the  inelastic  and  elastic-collision  contributions, 
sum  the  specific  energies  deposited  by  the  loss-cone  (PEI  and  PEE), 
ion-leak  (QEI  and  QEE),  and  CHEX  (CEI  and  CEE)  particles. 

e.(erg/g)  = PEI  + QEI  ^ CEI 
e^(erg/g)  = PEE  + QEE  ^ CEE  . 

3.  Compute  the  number  density  of  nuclei  in  N^,  and  O. 

n(nuclei/cm^)  = 2'|[N2j+  [02l(' " ’ 

4.  Compute  the  threshold-energy  parameter  used  in  modeling 
the  cascade  of  elastic  energy. 

E^  = -|39.  04^X2]  • 20.  48[02]  + 25.  40[Ol[/n  . 

In  this  equation  the  coefficient  of  [O  ] includes  the  replacement  of 
by  e^^^(O)  + e^,(0)  in  Eq.  (68)  of  HM-72b  (or  Eq.  (40)  of  HL-73b). 

5.  Compute  the  weighted  value  of  the  energy  removed  per 
collision  in  the  cascade  of  elastic  energv'. 

e(eV)  - lCjlN2|  - C2[02]  ^ C3[0|}  n 

where 
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Cj  = 2[9.76  + 2.37(1  - Fg)  + 12.  0] 

= 2[5.  12  + 1.  96  Fjg  - 12.0] 

Cg  - 1.96  F^g  + 12.0 
Fg  = 4/14 
= 5/14  . 

6.  Limit  the  specific  enerp^es  to  values  less  than  those  re- 
quired, according  to  the  model,  for  complete  ionization  and  (effective) 
dissociation  of  the  N^,  O^,  and  O densities  in  excess  of  1.0.  Define 

NHp(cnr^)  = [Ng]  ^ [O2]  + [O]  - 3.  0 


€. 

imax 


emax  1 HP 


where  the  conversion  factor  is 


. / ./-I  3 -1, 

e^(ergeV  cm  g ) = 


= 1.60x10  erg/eV 

The  fractions  of  destroyed  by  inelastic  and  by  elastic  processes 

are,  respectively, 

S-  = f-/e. 

1 1 innrx 

and 

R = ^ /f 
e e em;ix 
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H K.  + e > 1.0,  one  must  reduce  e.  and  e bv  the  ratio  l/(s.  - 

•^1  ®^e  1 e ■ 1 e . 

We  are  now  ready  to  deal  explicitly  with  the  inelastic  collisions. 

1.  Compute  the  electron  density: 

[e]  = e./(28.  5 . 


2.  Partition  the  charged  species  from  inelastic  collisions. 


IN^] 

[Og] 

IN"1 

[0^2 


1 

ZlNj] 

1.35 

n 

1 

2102 1 

1.  35 

n 

0.35 

21N2I 

1.  35 

n 

0.35 

2102I 

1.35 

n 

Lei 


Lei 


Lei 


Lei 


Lo'li  = ^Lel 

lo'  l - Lo‘l2  - Lo'l^  . 

3.  Record  the  corresponding  neutral-atom  densities  produced 
in  dissociative  ionization  in  inelastic  collisions: 

[N('^S)|  = [N^l 

10(^P)1  Lo‘)2  • 
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4.  Partition  the  enei's>'  from  inelastic  collisions  into  four 
modes  and  the  fraction  of  energy  in  each  mode. 


Electron  Thermal  Energy: 


E(eV/cm  ) = 1.  5 [ej 


'k  = S 


Dissociation  Energy: 


D.(eV/cm  ) - 9.  76  [N^  ] - 5.  12  [O^J2 


f'  = D.e. 
D 111 


Ionization  Energy: 


KeV/cm'^)  = 15.  58  [N"j  * 12.  06  [o"]  - 14.  53  [n"]  - 13.  62  [O^] 


Iladiation  Energy: 


'i  ' 'i 


Ij(erg/g)  = I 


R.(eV/cm  ) - I'  e./e,  . 

1 R 1 1 

We  now  deal  e.xplicitly  with  the  elastic  collisions: 

1.  Compute  the  total  number  density  of  particles  from  the 
cascade  of  elastic  energ\': 


Q 4 f e 
, -3x  KE  e 

V"''  > = 7,  i3e,  ■ 
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where 

E - c,  3 
^KE  " E^  ^ e 


2.  Partition  the  neutral  species  from  elastic  collisions.  The 
(fast)  particle  densities  from  the  elastic-energy  cascade  are 


21^2 1 

n 

n c 


2102I 


n 

c 


l°lc2  ■ 10 1 


cl 


The  (slow)  particle  densities  from  the  elastic-energ\'  cascade  are 

|N|„  = [N|^ 

10 1.2s  ' 10  1.2  ■ 

3.  Paiiition  among  excited  states  the  (total)  particle  densi- 
ties from  the  elastic-energy  cascade. 

IN(^S)|  = Fg[N|^  . iNl^g 
IN(2d)|  = (1  - Eg) 

I0(2p)|  = „ - F,„)  101^-101^2^ 
lO(^D)|  = Fjg[0|^,  . 
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4.  Partition  the  ener>^\'  from  elastic  collisions  into  four 


modes. 

Heavy -particle  thermal  eners>” 

K(eV/cm^)  = f e e 
KE  e 1 


K^(erg/g)  ^ K 

Dissociation  energy; 

D^(eV/cm^)  = 9.76[N]^  + 5. 12[0]^2 
Excitation  energy: 

X (eV/cm^)  = 2.37  (1  - F.)[N]  - 1.96  [Ol 
e y X 10  X 

X^j(erg/g)  = tj 


Radiation  energy: 

R(eV'Cm^)  = 12-j[N]^  + [0]^[  . 

We  can  now  combine  the  results  from  inelastic  and  elastic 
collisions. 

1.  Add  the  contributions  from  the  inelastic  and  elastic  colli- 
sions to  obtain  the  total  values  of  dissociation  and  radiation  energies. 

D(eV/cm^)  = D.  + D 
1 e 

Dj(erg/g)  = D 

R(eV/cm^  = R.  . R 
i e 
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Rj(erg/g)  - R • 


2.  Compute  the  total  fractional  partitions  of  enere^'  for  the 


six  modes. 


'k  = 'l  * ‘e> 


'd  ' 'l  °'<'i  ’ 'e* 


f^  = € X / €.  + e ) 
X lei  e 
e 


t,  = l/(e,  * £^) 


'r  ' 'l  "‘e'  • 


The  fractional  error  in  energy  consei'\’ation  is 


^ " ^X  " ^R  " ^ 

e 


3.  Compute  (and  collate)  the  changes  in  species  densities. 


6[N,1  = [e|.  n. 


-ile|  * .v| 


4 - 

6[n(^s)]  hn  h (1  ^ 


J 


2 2[N  ] 

6[N(^D)J  (1  - 

lot  - 2[o„] 

‘101=  --l^[e|.[0  |2  + 2-^n^ 

6[e]  = le] 

^ 1 2[N  ] 

^ I*-! 

. 1 2[Oj 

^^^2^  “ TT35  ~1T~ 

6[N^]  = 0.  35  6[N2] 

6LO”")  = 0.35  6[O2l  + i^he| 

4.  Compute  the  sums  of  changes  in  densities  of  nitrogen  nuclei, 
oxygen  nuclei,  and  charged  species  as  measures  of  consei'vation  errors. 

= 2 j6[N2J  f 6[N2 If  + 6[n(^S)|  + 6[N(^D)]  - 6[N^| 

X:^0  = 2 I6IO2]  . 6[02l}  + 6[OW  6[o"l 

- 6LN2I  t ^LOgl  + ^[N^l  4 6[0^]  - 6[el  . 

5.  Compute  the  sum  of  electron  thermal  and  0(^D)  excitation 
energies  per  newly-formed  electron. 

E^^^(eV/electron)  = (e  . 1.  96  |l02l  ^ 0.5lOlf^|/el  . 
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6.  Stoi’e  changes  in  cell  quantities  in  temporary  scratch  area 

(BUF2). 


BUF2(1)  = 6[e] 

BUF2(  7)  = 6[n"| 

BUF2(2)  = SLN("^S)] 

BUF2(  8)  = 6[0*  1 

BUF2(3)  = 6[n(^D)] 

BUF2(  9)  - E 

pe 

BUF2(4)  - 6[Oj 

BUF2(10)  = Kj 

BUF2(5)  = 6[N^] 

BUF2(11)  = 6[N2l 

BUF2(6)  = 6[02] 

BUF2(12)  = ^[O^] 
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8.  ORGANIZATION  OF  THE  HEAVY- PARTICLE 
RELATED  ROUTINES 

An  outline  of  the  portion  of  the  HAG  module  for  kinetic-yield 
partitioning,  heavy-particle  energ\'  deposition,  and  generation  of  conse- 
quent chemical  species  is  shown  in  Fig.  9.  The  interface  with  the  rest 
of  the  HAG  module  is  through  three  calls  by  Subroutine  PROMPG. 

The  first  call  (for  debris)  from  Subroutine  PROMPG  is  to  Sub- 
routine DEBRIS  which  seiwes  as  a control  routine  for  seven  subsequent 
steps.  Subroutine  HDPART  is  called  in  Step  1 for  the  partitioning  of  the 
kinetic  yield,  given  the  kinetic  yield,  event-point  magnetic  field,  and 
event-point  mass  density.  Step  1 is  accomplished  with  the  aid  of  Sub- 
routines LOSCON,  lONLEK,  and  CHXLOS.  A return  to  Subroutine  DEBRIS 
occurs  with  values  for  the  downward-moving  loss-cone  energy  (WLCD), 
the  downward-moving  ion- leak  energy  (WILD),  the  charge- exchange  en- 
ergy (WCHEX),  the  characteristic  magnetic  tube  radii  for  loss-cone 
particles  (R^^)  and  ion-leak  particles  (R-^)>  the  minimum  (v^)  and 
niivximum  (V2)  velocities  in  the  loss- cone  energy'  spectrum,  and  the 
residual  energy  available  for  UV  and  thermal  energies. 

If  WLCD  ^ 10  ^ MT,  then  Step  2 is  a call  to  Subroutine  CONSPC 
to  compute  the  histogram  energy'  distribution  for  the  loss-cone  particles 
and  Step  3 is  a call  to  Subroutine  BLINE  to  deposit  the  loss- cone  particle 
energies  in  a spatial  tube  centered  about  the  event-point  field  line.  Sub- 
routines CONJUG,  BFTELD,  and  INDEX  aid  in  the  field- line  tracing 
while  Subroutine  DEPO  does  the  parti(  le  energy  deposition  in  a grid  cell 
and  detects  the  stopping  locations  of  the  monoenergetic  energy  groups. 
Subroutine  XYZGEO  performs  the  necessary  transformation  from  tangent- 
plane  cartesian  coordinates  to  geocentric  cartesian  coordinates.  Wlien 
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Fiji,  9.  Simplified  Flow  Chart  of  Debris-Related  Routines. 
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control  returns  from  Subroutine  BLINE  to  Subroutine  DEBRIS,  every 
cell  that  is  intersected  by  the  loss- cone  magnetic  tube  has  a contribution 
to  its  specific  energy. 

If  WILD  ^ MT,  then  Step  4 is  a call  to  Subroutine  LEKSPC 
to  compute  the  histogram  energy  distribution  for  the  ion- leak  particles 
and  Step  5 is  a call  to  Subroutine  BLINE  to  deposit  the  ion- leak  particle 
energies  in  a spatial  tube  centered  about  the  event-point  field  line.  The 
operational  procedure  of  Subroutine  BLINE  is  identical  to  the  loss- cone 
case  except  that  the  heavy  particles  are  considered  to  have  mass  number 
16  rather  than  the  mass  number  27  of  the  loss- cone  particles.  After 
specific-energy  contributions  are  added  to  cells  intersected  by  the  ion- 
leak  tube,  control  is  returned  from  Subroutine  BLINE  to  Subroutine 
DEBRIS. 

If  WCHEX>  10"^  MT,  then  Step  6 is  a call  to  Subroutine  CHXSPC 
to  compute  the  histogram  energy  distribution  of  the  CHEX  particles  for 
later  use  by  Subroutine  CHXDEP.  Step  7 is  the  return  to  Subroutine 
PROMPG. 

The  second  call  from  Subroutine  PROMPG  to  the  debris  module 

-4 

is  to  Subroutine  CHXDEP,  provided  WCHEX>  10  MT:  otherwise  the 
call  is  bypassed.  Contributions  to  the  specific  energ>'  by  CHEX  particles 
are  assigned  to  every  grid  cell,  provided  the  areal  density  to  the  cell  is 
less  than  the  range  of  the  most  energetic  CHEX  particle. 

The  third  call  from  Subroutine  PROMPG  is  then  made.  Con- 
tributions to  cell  specific  energy  by  the  three  particle  sources  are 
summed  for  the  production  of  changes  in  the  chemical  species,  thermal 
energy,  and  excitation  energy.  These  latter  two  calls  from  Subroutine 
PROMPG  are  conveniently  made  just  prior  to  the  UV  photon  deposition 
since  the  necessary  cell-calling  logic  is  the  same. 
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Detailed  flow  charts  for  Subroutines  BLINE,  DEPO,  and  CHXDEP 
are  given  in  Figs.  10,  11,  and  12,  respectively. 


Fi{>;.  10.  F^low  Chart  for  Subroutine  BLINE(NXCEL,  RR). 
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Fig.  10.  (Continued). 
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Do  MMAX  points  on  circle 


DO  100  M = l.MNtAX  | 

II3 

Initiiilizv  cumulative  steppim*  distance  in  altitude  SDIIfcm),  flajts.  ^ 
heavy-particlc  ener^  bins  bffore  starting  t»  trace  neu-  field  line. 


FM  = M $ SDH  - 0.  0 
KDL  - 1 S DI.FACT  8. 0 
IFGNKW  = 1 $ IFLGV  - 0 
r-DO  30  L - l.IWK 
I KL(L)=^F.KfL) 

30  CONTl.VLF 


Start  non  event-point  field  line.  Determine  position  in  tancent-plane 
coordinate  system  (X.  Y.  V)  and  transform  to  ge‘>cen*nc  spherical 
coordinate  system.  Determine  equatorial  altitude  HUM  for  field  line 
through  Point  M. 


PSIM  = 2.  MM-FM  SWIF2  - SHIFT 
X = SMALI.R  • COS(PSIM)  $ Y SMALLR  • SI.S(PSIM) 

Z 0.  S HTPO  - 1.E7 

CALL  XYZGEOlHTPO,  GCTPO,  GLTPO,  X,  Y,  Z,  HAM.  GCM,  GI*M) 
HAM  KM  ^ H.AM  • l.K-5 

CALL  BFIFLI)(HALFP1-GCM,  GLM.  HAMKM.  lU'AL,  DIPANG. 
DKCANG.CO.SCHI) 


GROM  ^ (RE  - H.XMV'U.  - COSCVH'*2) 
HOM  = GROM  - RE 


Set  starting  values  of  field-line  tracing  parameters  to  those  at  event 
point. 


50  cev  = OCX  s GLV  = GLX  S HAV  = HAX 

GRV  - RE  ♦ HAV  $ BV  PP(-l)  $ DIPX  ==  PP(51 
DF.CX  = PP(6)  $ CSCHIX  - PP(7) 


Determine  initi.il  altitude  step  DM.  to  be  taken  Liter  and  held 
con.stant  during  traversal  of  fir.*:t  cell.  Provide  for  non- zero  step  if 
p< tint  13  at  c,'«|uator.  

Dll  DL  • Sir,’(  DII>X  1 
IFtDH.EQ.  0.0>  DM  = 1.K4 

t 

(.0 

Fig.  10.  (Continued). 
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Fip.  10.  (Continued). 
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-52-!- 


If  ceU-travcr^al  Conipl»*»ed.  »*  • t “0  1 1 con  |ju?f  pirh  d*  ;'  Mt  ••tnT^’v,  aiui 

t*ithor  renmn.jtp  track  if  appr  •pruio  * r pMaMi-t.  r » s f •!  nt-Kt-n-n 

tra%'ers:il.  II  cell-traversal  m-t  c<jinpU'»ed.  continue  stpppir.i*  t v e nnp  to  b'j. 


L-  - . 

L 


Fig.  10.  (Coiitinuod). 
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Compute  field- line  leriidh  from  entrance  to  exit  of  cell,  «|ua'onul  value  of  ft-  lj  I'.hy, 
equatorial  radius  RKQ,  and  areal  mass  densitv  ZKTA  alone  livid  line  from  lapurtixi- 
mate)  en>  ranee  to  (approximate)  exit.  Call  Sutiroutine  l>Kl*t)  to  conipute  envi'!  v 
deposited  in  cell  NSCt'I,.  ff  cell  d is  outs.de  of  (tnd.  avoid  deieiMtino  enerv-,-  in 
nearest-cell  iii  prid  by  pomp  to  75.  Comjiute  piocontric  i.idiu*  (illl'fiOT  to  h.  itoiu  of  < 
cell  NSCKL,  column  number  NtSCOl.  for  cell  NrCKL,  steradian.il  volume  SHVULS  of  j 

cell  NSCEI.,  factor  F.PS'.’.AS  to  convert  from  eV  cell  to  erp'c.  and  enerpy  deposited  ; 

in  cell  Nb'CFL  due  to  elastic  and  inelastic  collisions  (both  on;  p and  r-V.  cell) 


70  BF.Q  - BKBl  • BS 

RF;Q=  (l.EelS  • .MUO/BEQ)’ *7111111) 

ZETA  = S12  • IlEQ  * RliOKKiN'SCEL) 

CALI,  DEPOlZEIA.IJIV.ANEE.A.N'EI) 

IF  (IFLCS.EQ.  1)  «)  TO  75 
CUSBOr  * GRINSCEI.  - 1) 

I,K  = MODfNSCEL.N’CELZ) 

IF  (LK.EQ.  1)  GRiBOT  - RE  » HABOT 

NSCOL  - (NSCEL  - 11,  NCELZ  » 1 

SRV01.S  = (GR(NSCEI.1**3  - GRSBOT *3)*THIRD 

EPSMAS  = EPSl.N  (SRVOLS  * STER(NaCOL)  • RIIOKK(N.SCEL)) 

QEEINSCEU  = QEEIN.'tCEL)  - ANEE'EPSMAS 

QEI(NSCEL)  =QE1(NSCEL)  - ANEl'EPSMAS 

ENEE(NSCEL)  = EN’EE(NSCEL)  » ANEE 

ENEl(NSCEL)  = ENEI(NSCEL)  * ANEI 


100. 


Fiff.  10,  (Continued). 
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ITsf  r:tr.''t'-c[i*'r^r\-  parano-ti  i •'! 
oxyp,ea  in  niti  t»F,en  lor  ».>n-U-  li. 
p.mi<  h 

^ ! ^ 

E lOAi)  0 Vi  ► 

/ 2 

L'se  1 .in;:i.*-(,i.ci  I'V  p.iranu-tera  (»f  ' 
.ilutmmim  m lotiort  n for  lus.s-  > 

cone  p.ulirles  j 

ALPHA  Al.PilAI  $ IlKTA  HU  A 1 

AM’IIA  Al.l'ilA25  m TA  ^ llhTA’' 

('  -inputt'  rediic»*J  r.itit’t*  ZLTAK  uilliin  cell  NSCKL 

I<A  liKIA  AM'itA 
FAC  iUI^  0.  r>  ' m.  lA  • hA  .'*M 
AM  1.  l.\  > Al<2  An  • Ali 
/FTAV  - V’lllA  * FAClOK 


Compute  cell-cxit  energies  (e\')  of  particKs  :n  |u.<s-c**ne  i»r  i-m-loak  enonrv  hiJ'topram,  thfjr  traversed  * 
rant:o,  and  — if  nut  exitnif'.  — tin-ir  .•'tnppit,.-.  altitudr-s  and  C'  ■ '.-.i .iphic  coordinates.  j 


Compute  reduced  r:in<;e  7KTAf  1 c-f  incornnut  p.'irtfclc  with  rcductd  enerjzv  KlKi-.I>  Initialize  exit  energvi 
Efc'(l)  of  Group-1  pvirticle.  { 


DO  20  I l.IV.K 
HRKD  R\  • SQRT(KI.'I)^ 

ZKIAH  - KUiU>  - AIa^u.1-0  • r.lRKD) 
t£(l)  - 0.0 


t£(I}  - 0.0  J 

^ — : : 

I Compare  particle  ran^o  ZKi.AKI  wilii  cell  path-Undh  7K1‘.\F 


Compute  cell-exit  enerjt\-  of  U 
particle.  j 

ZETFlll  ^ 7KTAF  EillFD  i 

TERM  = SQia’<2*ZE7Fin-l  ' - 1 ; 
FEFOTK  2n*FlRFD)iZi:ii:iK  | 
-TEftM)  L'lRKD  ' 

^ ' 

KKKOTE>0.  001  ^ 

- -1 

EF(I)  = AH2[a:mFD)2  - 2’E;RM>' 

>(l.FIlU-;D)^lFK\tJ  i 


DFI.ZF  l 7FTAFI  - /FTAl- 
FilQm  I-/F1/OKL7.KT.IFG) 
FK(I)  AH2  - .rSQll'l  )^ 


Zr  iAF.  GE.  ZKIAI  I 


t 

DZFTA  ^ ZKIA  I 


DZHIAa) 


= 0.0  ;Ui  J 


El.(l)  - 0.0 

“if  |f 


Compute  altitude  and  p.ef'p.r.iphic  c«x»rdi- 
n.tles  of  stoppir.j,  pujr.t  of  particle.  Com- 
pute ranpe  ZFT.Al  • f inromir.!:  particle 
with  eneri5\'  III.  I).  F stimat'-  decrement 
in  attitude  DKLIIl-wM  Fy  sc.ilinc  with  ratio 
of  r.inge.s.  SuFti  act  il«’cr»  nient  from 
cell-entrance  altitude.  Call  CUNMl/G  t ' 
ol>tatn  4;eo{;raphic  ivrth  l.itUude  and  o.tst 
I 'nr.itude  of  stoppin!:  j^iint.  Convert 
Hcui:iaphic  col.ttduil*-.  C >nvr  rt  altitude 
from  km  to  cm. 


ZKIAI  ZKlAll  lACrOR 
DZI  I A n ZIIAJ 

DLl.HKM  (HA.-KM  - MAVKM ‘'ZETA 1 
7.W  1 A 

! ii\>;km  ha.skv  - nEi.nir^! 

I CALI.  (’ON.U  G'fK'L-.Gl  ; . HASfrvi. 

I HVI'KM,  I.0.C.CU;,GLE.S12. 

I DEnn 

; GCi:  - HAl.m  - GCLI. 

I nvr  iivEKM  *11^ 

! CAKDIV  (.'NKP’PIV 

I WlUlMC.ll  . I.QMvl)IV.GCF. 

I GH..H\i; 

, U l-ORMAT 
I C/)  T(.)  20 


20  CONIINIT 


FiR.  11.  Flow  Chart  for  Subroutine  DEPO(ZETA.  DIV,  ANEE,  ANP:i). 
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RtnURN  ^ 


Fip.  11.  (Continued) 


Compute  cell-entrance  and  cell-exit  enertnes  (eV)  of  particles  in 
CHEX  enerRV'  histujtram  and  their  traversed  ranges  tg  cm**).  Com- 
pute reduced  range  ZETA.Fl  of  particle  vuth  reduced  enerpv  EIRED. 
Initialize  cell-entrance  and  cell-exit  energies  EKl)  and  EOiI). 


?0  CONTINUE  I 


Flow  Chart  for  Subroutine  CHXDEP(NCEL). 
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Fig.  12.  (Continued^. 
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9.  UV  DEPOSITION 


The  total  computation  of  UV  deposition  requires  Subroutine 
EUXFIT  for  finding  the  energy  in  each  of  the  five  Group- U subgroups, 
Subroutine  PLINE  for  tracing  ray  paths  from  the  event  point  through  the 
center  of  each  of  the  cells  in  the  HAG,  Subroutine  PINT  to  compute 
integrals  of  mass  density  and  species  number  densities  along  these  ray 
paths.  Subroutine  HEDGE  to  find  the  radius  (B-edge)  within  which  only 
atomic  ions  are  formed  and  to  compute  line  integrals  of  species  densi- 
ties within  the  B-edge,  Subroutine  PHEAT  to  compute  the  photon  fluences 
of  the  Group- U subgroups  and  to  compute  the  total  pressure  due  to  the 
Group- U energy  deposition  (and  to  combine  it  with  the  original  pressure 
and  that  due  to  Groups  X and  H),  and  Subroutine  PCHEM  to  compute  the 
prompt  chemistry  associated  with  the  Group-U  deposition  (and  to  com- 
bine these  results  with  those  due  to  Groups  X and  H). 

9.  1 COMPUTATION  OF  RAY  PATHS  THROUGH  THE  HAG 
(SUBROUTINE  PLINE) 

9.  1.  1 Int  roduction 

Subroutine  PLINE  traces  ray  paths  through  the  HAG.  A ray 
path  starts  at  the  event  point,  passes  through  the  center  of  a selected 
HAG  cell  regarded  as  a target  cell,  and  ends  as  the  ray  leaves  the  target 
cell.  We  must  compute  the  length  of  each  of  the  line  segments  formed 
by  the  intersection  of  the  ray  with  the  inteiwening  cell  boundary  surfaces. 
I Target  cells  are  selected  by  Subroutine  PROMPG  bv  systematically 

looping  over  all  the  cells  in  the  HAG. 
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9.1.2  Transformation  from  the  Quadrupole  Coordinate  System  to 
the  Associated  Geocentric  Cartesian  System 

In  Volume  16,  we  have  intrfxluced  the  quadrupole  coordinate 
system  and  given  many  of  the  details.  Here,  as  shown  in  Fig.  13,  we 
need  to  convert  from  quadrupole  coordinates  (quadrupole  colatitude  o 
and  quadrupole  east  longitude  6)  and  geocentric  radius  r to  the  corre- 
sponding geocentric  cartesian  system  in  which  the  Z axis  is  through  the 
quadrupole  north  pole  and  the  X axis  is  through  the  quadrupole  west 
pole. 

The  required  transformation  equations  are 
.X  = r sin  a cos  9 
y = r sin  a sin  6 
z = r cos  a , 

where  sin  a is  obtained  by  first  applying  the  law  of  sines  to  triangle 
P-QN-QW  to  obtain 

sin  a - sin  o sin  90  /sin  A 

./  2 \* 

= sin  o/ \\  - cos  A j 

and  then  the  law  of  cosines  for  angles  to  obtain 

cos  A = - cos  o cos  6 


so  that 


sin  a = sincp  ll  - cos  o cos 


2 2„\‘ 
II  - cos  O cos  9 1 


We  also  have 


cos  a = SIGN(1.  0,  cos 


1 - sin  a 


I 
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9.1.3  Cell  Shape 


i 

i 

I 


tr 


!j 


i? 
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As  explained  in  Volume  16,  for  a cell  in  the  HAG,  the  top 
and  bottom  surfaces  are  geocentric  spherical  surfaces  and  the  four 
sides  are  meridian  planes. 

9.1.4  Basic  Formulas  for  Slant  Range  from  a Given  Point 
to  a Cell  Boundary 

There  are  three  basic  formulas,  each  one  corresponding  to  the 
intersected  boundary  being  (a)  a top  or  bottom  boundai'y,  (b)  a quadru- 
pole  colatitude  boundary,  or  (c)  a quadnipole  east  longitude  boundary. 
The  first  case  affords  more  complication  than  either  of  the  remaining 
two,  each  of  which  is  similar  to  the  other. 


The  quantities  we  have  as  input  are: 

a.  Event-point  quadrupole  coordinates  and,  hence,  cartesian 
coordinates 

b.  Target-point  quadrupole  coordinates  and,  hence,  cartesian 
coordinates 

c.  Lateral  cell-boundary  quadrupole  coordinates  and  top  and 
bottom  geocentric  radii. 


In  each  of  the  three  types  of  boundary  cases,  the  ray  from  the 
starting  point  (usually  the  event  point,  denoted  by  r^)  through  the  target- 
cell center  (denoted  here  by  r^)  is  specified  by  its  direction  cosines. 


cos  a = (Xo  - 

X , ) r „ „ 

3 

1 13 

cos  i3  = {>-3  - 

cos  y = (Z3  - 

where 


13 


(.X3  - . (V3  - y^)^  1 (Z3  - Zj)2 


77 


I 


9. 1.4.1  Slant  Range  of  a Ray  Intei'secting  a Cell  Top-  or 


Bottom  - Boundary . 


Consider  a Point  2 with  geocentric  radius  on 


the  ray  shown  in  Fig.  14.  The  angle  <5  between  the  yectors  r^ 

— > ■Hi 

and  r^2  ~ ^'2  ~ specified  in  terms  of  the  direction  cosines  of  the 

unit  yectors  f ^ and  f the  relation 


^1  ^1  ^1 

r.  • r._  = cos  6 = — cos  oc  - — cos  ^ + — cos  y 

1 12  r^  r^  r j 


Xj  cos  a + cos  j3  -r  cos  y 


2 

>'l 


From  the  right  triangle  OPP^  in  Fig.  14,  we  haye 


(r^  sin  6)  + (r^  cos  6 - r 


1 


12' 


r 


2 

2 


DSRR  H r^2  " - - FGR  WS2 

where 

WS2  = (S2)^  + (GR2)^  - GR12 
S2  = r^  cos  6 
GR2  = i'2 
GR12  = I’j 


and  FGR  is  chosen  to  select  the  positive  root  or  the  smaller  of  the  two 
positive  roots. 
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9. 1.  4.  2 Slant  ftange  of  a R^iv  Intersecting  a Cell  North-  or 
South- Boundary  (Specified  by  Quadrupole  Colatitude).  In  this  case  we 
are  given  the  quadrupole  colatitude  (^2)  of  the  Point  2.  In  Fig.  15,  the 
line  segment  P^P2  projected  onto  the  YZ  plane.  Thus, 


sincp„  y 

tan  ^ — 

2 cos 


yi  - ■'12  ^ 

z,  - r,2  cosy 


or 


DSRC  = r^2 


Vj  cos  <^2  " sin 
cos  y sin  <^2  - eos  cos 


9.  1.4.3  Slant  Range  of  a liav  Intersecting  a Cell  East-  or 
West- Boundary  (Specified  by  Quadrupole  East  Longitude).  In  this  case 
we  are  given  the  quadrupole  east  longitude  (O2)  of  the  Point  2.  In 
Fig.  16,  the  line  segment  P|P2  is  projected  onto  the  XY  plane.  Thus, 


sin  60 

tan  6 = — 

2 cos  9, 


>^1 


r^2  00s  (i 
r^2  ^os  0; 


or 


DSRL  r^2 


Vj  cos  ^2 
cos  a sin  92 


- Xj  sin  92 

- cos  /i  cos  «2 


9.  1 . 5 A Caution  for  Energy  Deposition  in  the  Event- Point  Cell 

The  nominal  cell  heights  for  both  SAIHYD  and  NRLHYD,  under 
' 't!  initial  and  rezone  conditions,  are  given  in  Tables  2 and  4 of  Volume 
! . 1 r -m  these  tables  it  is  seen  that  the  cell  thicknesses  for 

t 't  the  cells  (and  especiallv  for  the  lower-altitude  cells)  are  much 
ir.iih  r than  the  typical  horizontal  dimensions  of  the  columns  wliich  mav 
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cos  a 


Fig.  16.  Geometry  Used  in  Deriving  Slant  FUinge  of  a R;iv  Intersecting 
a Cell  East-  or  West-Boundary  (Specified  iiv  Quadrupole 
East  Longitude  62)-  The  segnient  P1P2  >s  projected  onto 
the  XY  plane. 


be  on  the  order  of  100  km.  For  a small-energ\’  event  in  the  lower- 
altitude  cells,  it  would  be  quite  possible  to  put  more  energ>’  into  a cell 
than  was  available  from  the  event,  owing  to  using  a path  length  that  was 
unrepresentatively  small.  We  tried  an  interim  procedure  to  relieve  this 
situation,  but  the  undesii'able  side  effects  led  us  to  delete  the  procedure. 
Thus,  we  alert  the  user  to  this  potential  problem  which  deserves  review 
and  more  study. 

9.1.6  Directional  Flags  and  Indices  of  Next  Set  of  Boundaries 
to  be  Intersected 

In  Subroutine  PLINE  we  use  the  following  directional  flags: 


Index 

Value 

Play  Moves 

IINC 

t 1 

southward 

• -1 

northward 

.TING 

) 

eastward 

* -1 

westward 

KING 

) 1 

upward 

t -1 

downward 

The  initial  value  of  KING  is  saved  as  the  variable  KINCO. 

The  indices  I,  J,  K are  used  to  denote  the  ne.xt  possible  set  of 
boundaries  to  be  intersected  by  the  ray.  In  starting  from  the  event  cell 
(lEC,  JEC,  KEC),  these  boundary  indices  are  imtiallv  defined  as 

I = lEC  - (1  - IfNC)  2 
J = JEC  - (1  - JINO/2 
K = KF,C  - (1  - KINC)/2  . 

A value  of  1 (instead  of  0)  for  the  flag  [END  denotes  that  the 
last  step  in  the  line  integral  has  been  reached. 


The  index  II  is  used  to  denote  the  number  of  boundaries  inter- 
sected by  the  ray  from  the  event  point  to  the  exit  boundary  of  the  target 
cell- 

9.1.7  Flow  Charts  of  Subroutine  PLINE 

Figure  17  is  a simplified  flow  chart  of  Subroutine  PLINE  which 
includes  calls  to  Subroutines  PINT  and  BEDGE.  Figure  18  is  another 
simplified  flow  chart  of  Subroutine  PLINE  that  gives  information  only 
about  the  computation  of  the  ray- path  segments. 

9.1.8  Miscellaneous  Functions  of  Subi^outine  PLINE 

Other  miscellaneous  related  tasks  performed  by  Subroutine 
PLINE  include  the  following: 

a.  Set  Group-U  subgroup  index  (IGROUP)  to  5. 

b.  Zero  integrals  (SNR2(L),  FLUX(L),  L=1 , 5)  computed  in 
Subroutine  BEDGE. 

c.  Zero  integrals  (SN2INT,  S02INT,  SNTNT,  SOINT,  RHOINT, 
02INT)  computed  in  Subroutine  PINT. 

d.  Set  flag  (lUFLAG)  denoting  whether  or  not  the  B-edge  has 
been  passed. 

e.  Zero  integral  (SR)  for  distance  from  event  point  to  cell 
entry-point. 

9.  2 COMPUTATION  OF  INTEGRALS  OF  SPECIES  NUMBER 

DENSITIES  AND  MASS  DENSITIES  ALONG  THE  lUY  PATH 
(SUBROUTINE  PINT) 

For  the  line  path  SRN  calculated  by  Subroutine  PLINE,  Sub- 
routine PINT  integrates  the  species  number  densities  and  mass  densities 
as  follows; 
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Determine  direction  cosines  of  niy  from  strut-i'  i.nt  throucli  tarei  t rwint.  j 


Set  directional  flass  iitdicatinR  nhetlier  ray  poes  sooth  or  north,  east  ..r  ,-  .-1.  and  op  ; 
or  down.  Set  imtice.s  of  next  boundarie.s  to  lie  intersected. 


Compute  first  Set  of  slant  ranees  IDSHC.DSRL,  and  D.--nK)  from  even!  p...i,;  to  t-v.  -. 
cell  boundaries  (c.ilatitude,  loiipituiie.  and  geocentric  radius).  ; 


130 


End  of  line  into^ral 


no,  270 


DSRR  sh  ortest 


r, 


19of 

Sfc;  ==  DSSR  ^ 
CALL  IMNT 

Compute  DSRR  [ 
to  Opposite  [-« 
boundary  of  cell  i 
J 


DSRt,  shortest  ►j^DERC  shoi  tr.M  ] 

t -r 

2CoT 

f 

[ ^S!L\  = DSRI.  1 

1 SR.S  - DSKC 

CALL  PLNT  1 

1 CALL  l'l\T 

j 

i 

Compete  DSRL  ] 

( — 1 

1 Compute  DSKC  \ 

to  opposite 

■ to  opposite  j 

' 

boundary  of  cell ; 

— i 

j boundary  of  cell  , 

1 

1 

330 

Compute  vertical  index  K of  cell  in  nrw 

column  traversed  In 

ray.  ; 

390^ 

l^'ALI.  fiKDGE  j 


380  Y 

I’a.st  H-edge 

].j.— 


t 

End  of  path 


•110 


I'ln  ■ 


Fifj.  17.  Simplified  Flow  Chart  of  Subroutine  PLINE,  Includinii 
Calls  to  Subroutines  PINT  and  HEDGE. 
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Compute  direction  cosines  (CSALP,  CSBET.  Jnd  CSGAM)  of  ray  fi  viin  event-point 
(QLX,  QCX,  CriX)  throufth  taryet-cell  (IP.  -U’,  fiP)  center  and  S2  and  .S4. 


Set  directional  flacs  (IC.'C.  JINC.  KINO  and  indices  (1,  J,  K)  of  next  txiundanes  to  be 
intersected. 


< QL(.IP)  > QLX 


J - JEC  - (I  - .IINC)  2 i 
IfNC  ^ 1 I 


I = EC  - U - IINC.'2) 
KINC  = 1 


^ 

KEC  - U - KINC) '2 


i JINC 

= -I  1 

-]  IfNC 

I 

J 

■|  KINC 

■ J 

Save  initial  value  of  KING.  i 

KINCO  --  KING  I 

Compute  WS.  WSl.  WS2.  j 


Fig.  18.  Simplified  P’low  Chart  of  Subroutine  PLINF.  Onlv  .'iteps 
pertaining  to  ray-pnth  segments  are  shown. 
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Compute  first  set  of  slant  lanttes  from  eveiii-point  to  event-p*Hnl-cell  lx*i 
their  extcnsums). 


Compute  DSHC 
DSUC  > 0 


DSRC  1.  K20 


T 1 

. _ :j 

...J 

r"* 

Comput 

0 DSKI,  j 

<"  DSUL>0 

9f.T 

[ FRG  - - 1 • Ir-'r  f <•  0) 

S2.0  \ 

:.t-7 

I Compute  \V>2  j 

V - - 

WS2  2 0 7 

ioo|'_ 

KIN'C  - 1 

‘-’I.-,: 

iiof  __ 

[ WS2  = Max(0,  WS2) 


>E-  Dsri 


DSriL  1.E20 


['UG  - I -•  (COS  6 > 0) 


KIKC  = 1 , K = K • KINC  j 
FUG  - 1.  Compute  WS2  i 
__ 


T>20  _ 

^ DSrtll  = - E2  . KKG  ^^82  J 


llSItK  - M;ix(lKHK,  1) 


Fiji.  18.  (Continued). 
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SKN'  = Mjr.lD.'-RC.nSUI, 
DSURj 


K«'turn  to 
PHOMPG 


Compute  suhs<-«|ueiU  sets  of  sLint  r.int;cs  tu  iniervomr 
nes  (or  their  exunsiuns). 
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SN2INT 


S02INT 


SNINT 


SOINT 


RHOrNT 


02INT 


1^2 1 


[Oj] 


dr  , 


dr  , 


-2 

cm 


-2 

cm 


-2 

cm 


-2 

cm 


s/cm 


2 


g/ cm 


2 


The  variables  for  these  integrals  are  zeroed  in  Subroutine  PLINE.  Over 
a ray  segment  between  SR  and  SRN,  the  corresponding  contributions  to 
the  integrals  are,  respectively,  DSN2I,  DS021,  DSNI,  DSOI,  DHROl, 
and  D02INT.  SPECIE (1)  denotes  the  concentration  of  N^,  O2.  N.  and 
O,  respectively,  for  I = 1,2,3,  and  4. 

When  the  properties  of  a paidicular  cell  are  needed  to  perform 
the  integrals  Subroutine  PINT  calls  Subroutine  ECRD  and  transfers 
quantities  2 through  11  from  large  core  memory  to  small  core  memory. 
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9.  3 COMPUTATION  OF  THE  ATOMIC- ION  REGION  (B-EDGE) 

FOR  GROUP- U PHOTONS  (SUBROUTINE  BEDGE) 

9.3.1  Introduction 

As  Subroutine  PLINE  computes  the  ray  path  from  the  event 
point  to  a target  point,  Subroutine  BEDGE  is  employed  to  determine  the 
radius  (SRR(2))  of  the  edge  (B-edge)  of  the  region  within  which  Group- U 
ionizing  photons  are  very  stronglv  absorbed  and  produce  entirely  atomic 
ions.  Subroutine  BEDGE  also  computes  (a)  the  line  integrals,  from  the 
event  point  to  the  B-edge,  of  the  02-mass  density  (02U)  and  of  the 
species  densities  of  N2  (SN2U),  O2  (S02U),  N (SNU),  and  O (SOU),  and 
(b)  the  number  of  ionizing-photon  absorptions  per  steradian  that  are  re- 
quired to  produce  atomic  ions  alone  for  Group-U  Subgroup- L (SNR2(D). 

9.3.2  Simplified  Description  of  the  E ssence  of  Subroutine  BEDGE 

To  simplify  the  description  of  the  essence  of  Subroutine  BEDGE, 
we  introduce  some  terms,  not  all  of  which  are  used  in  the  routine. 

PHOTON(L)  = Number  of  Subgroup- L photons  emitted  per  steradian 
from  the  source. 

PHOTAB(L)  = Number  of  Subgroup-L  photons  per  steradian  used  in  de- 
veloping the  radius  of  the  B-edge  up  to  its  current 

radius  R . 

c 

ABSORB  = Number  of  ionizing-photon  absorptions  per  steradian  re- 
quired to  increase  the  radius  of  the  B-edge  from  its 

current  value  R to  R. 

c 

2 

XMUBAR(L)  = Mass  attenuation  coefficient  (cm  g)  for  Subgroup-L 
photons  in  Cell  NCELS  with  neutral-particle  mass- 
density  RHON(NCELS). 

4 

X]  SPECIE (1)  X SIGUIL,  I) 

I 1 

~RHON(NCELS) 

(I  1,  2,  3,  4 denotes  N2.  O2'  Ni  ^^) 
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EVEN  = Number  of  absorptions  per  cubic  centimeter  required  to 
get  atomic  ions. 


2 4 

• SPECIE  (I) 
1=1  1=3 


= SPECIE(I) 


EVNPG  = Number  of  absorptions  per  gram  I'equired  to  get  atomic 
ions. 

= EVEN/RHON(NCELS) 


If  we  want  the  Subgroup- L photons  to  produce  all  atomic  ions 
at  radius  R,  then  the  condition  to  be  satisfied  is 


PHOTON(L)  - PHOTAB(L)  - ADSORP 


X XMUBAR(L)  = EVNPG  . 


2 

Now,  multiply  both  sides  of  the  equation  by  R /(XMUBAR(L)  x EVEN) 
and  define  a new  function  EDGE(L)  to  be  the  difference  between  the  right- 
and  left-hand  sides; 


EDGE(L)  = 


SmObA^T  ■ ^PHOTON(L)  - PHOTAB(L)  - ABSORB) 

EVEN 


A positive  value  for  EDGE(L)  is  the  steradianal  volume  in  e.xcess  of  that 
for  whidi  the  Subgroup- L photons  can  produce  all  atomic  ions.  A nega- 
tive value  for  EDGE(L)  is  the  steradianal  volume  beyond  radius  R that 
the  Subgroup- L photons  can  produce  all  atomic  ions.  A zero  value  for 
EDGE(L)  means  that  the  Subgroup-L  photons  have  produced  atoniic  ions 
exactly  to  the  radius  R. 

If  we  recognize  the  fact  that 


ABSORB 

EVEN 


I 
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then  we  can  write  EDGE(L)  as 


EDGE(L)  = A„R^  t A_R^  - A 
3 2 o 


where 


A3  = 1/3 

Ag  = l/[RHON(NCELS)  x XMUBAR(L)] 

r3 

, _ PHOTAB(L)  - PHOTON(L)  c 

o EVEN  ■ 3 

The  equation 

3 X EDGE(L)  = 0 


is  solved  for  the  radius  R by  using  the  standard  Newton- Raphson  itera- 


tion formula: 

EDGE(L,  R.) 

R = R. 5 

1 d[EDGE(L,  R^]/dR. 


9.3.3  Flow  Charts  of  Subroutine  BEDGE 

Figures  19  and  20  are  simplified  and  detailed  flow  charts  of 
Subroutine  BEDGE. 

9.  4 COMPUTATION  OF  GROUP-U  SUBGROUP  FLUENCES  AND 
PROMPT  HEATING  (SUBROUTINE  PHEAT) 

For  each  cell  in  the  HAG,  Subroutine  PHEAT  computes  the 
photon  fluences  of  the  Grnup-U  subgroups  and  the  total  pressure  due  to 
the  Group-U  energ>'  deposition  and  then  combines  this  pressure  with  tlie 
pre-event  pressure  and  the  pressures  due  to  Group-X  and  -H  energy 
depositions. 
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^/iCRdi  ^^•T  IVLSOLU^I' 'N  M’’  'MAh’ 
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Compute  EUGI-.{S11N,  L=5) 


SRN  > S1U(L^5) 


Estimate  D-edi’c  radius  (SKI) 


Newtoa  iteration  method  for  f<-edno  radius  (SKI(L)) 


SRK<  SiU  • i>RS 


J Kisection  method! 
1 for  nev.  SRI 


(sra  - SRI)  SRI  0.001 


> 20  iteiutif'iis 


SRI  Min(SKl.SRl) 


SRN  > SRI  11) 


Fig.  19.  Simplified  Flow  Chart  of  Subroutine  BEDGE. 
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Initialize  absorptions  (sr~^)  up  to  radius  SUB  for  Suiju«'>np  L | 

3 I 

RHOBOT  = SNU2(L)  $ SRB3  - (SUB)  $ Uclour  if  sou rce-em-rr.y  1.  not  positive,  j 


WU(L)>  0 


M =■  0 S SRI  = SRli  j 


Save  EDCiK,  SUM,  and  SUN.  Compute  new  KDGK  at  radius  SRN.  Save  p:DGE0,  SUM, 
and  SRN  for  use  in  Bisection  Mflh<»d,  if  necessary. 


RHOTOT  = UHOMfJT  ^ DSNU2  S faJGEO  = EDGE 

= kO;n 


eugUTl] 


FO  = EDGED 
XA  = SRB 
XB  = SRN 


passed  B-edt;e 


10  T Passed  M-edge 

p g.  i.  

Set  coefficients  (AO  and  A2^  of  U and  R in  cubic  f 'r  iterative  solution  of  equation  \ 
EDGE  - 0.  Save  Group-E  LDGP'  for  radius  SitS.  Estirmtc  Group-L  B-edge  radius 
(SRI)  by  linearlv  interp  ilaiint:  between  SRI5  (corresponding  to  EDGK-value  of  EDGED)  ' 
and  SRN  (corresponding  to  EDGE-valuo  of  EDGE).  Zero  iteration  ci*untcr  (M).  1 


AO  * (3  'EVEN)  [KHOBOT  - WUlL).  KHGUiL)]  - (SHB)*^  S A2  * 3 'iRMOSEM  v XMffUH)  ' 
BEN(L)  = EDGE  S SRI  = (SRB  x EDGE  - SRN  • EDGED)  (EDGE  - EDGED)  $ M 0 j 


Use  Newton’s  iteration  method  to  .‘^olve  EDGE  equation  for  next  estin^te  of  B-cdge 
radius  (SKI).  = derivative  of  EDGE  with  ro.''p*T‘  t.,  radius. 

WS  = (3  X SKI  * 2 ^ A2)  ' SRI 

U S D • - — - 

n ''  _ 

F^valuate  Newton  formula.  If  SRI  n>'i  b(‘tween  SUM  and  SK.N,  use  Mi'^ertion  Method  to 
estimate  SRI.  SiU  SRI  - [(SUli-  (SRl  . A2  • AOj  ws 


^ H ' 

70V^  / 


SUf  < SRB  or  SKJ  > SUN 


Increment  counter.  Test  for  C'*n-  H 
vergence.  Eimit  iterations  to  20.  l 
Must  use  snuiUer  of  last  two 
iterations.  .M  .\1  * 1. 


I SRI  - SRI  SRI  < I.  E-3 


fbsci-tiun  Mrthi>d  | 

Compute  new  EDGE,  called  F2.  ' 
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UHOSUl  lUiOMUT 

. FVF.NliSmi^  - SRH3)  3 


3 

fiSltn^  FVNPG 

p;VEN 

[ XMEBAK 

M > 20  ©■>'  

- YJ5 
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SUl  « SRI  I- \ 
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RHOsm  I 
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Fig.  20.  (Continued). 


XA  - SUt  $ PO 


SRI  0.  S(XA  ♦ XB) 
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Fig.  20.  (Continued). 
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9.4.1  Group- U Energy  Deposition  and  Fluences 


In  computing  the  Group-U  energy  deposition  (erg/g)  in  a target 

2^ 


cell  and  the  Group-U  photon  fluences  (photons/cm  ),  we  treat  three 
cases: 


Case  1 (B-edge  radius  (SRR(2))  smaller  than  target- ceil  entry- radius 
(SR)) 


We  use  the  difference  of  the  incident  and  exit  values  of  the 
steradianal  energy  fluence  for  each  Group-U  subgroup,  divided  by  the 
corresponding  steradianal  mass  for  the  cell.  We  treat  dissociating 
(Subgroup  1)  and  ionizing  (Subgroups  2-5)  separately. 


DAIX(l)  = 


For  Subgroup  1, 
WU(1) 


-XMU1X02IN 


(p/3)  (R^  ^ - R^  ) 
^ out  in 


- e 


•XMU1X020UT 


where  XMUl(cm  /g)  is  the  mass  attenuation  coefficient  of  for 
Subgroup-1  photons,  02IN  and  020UT  are  the  areal  mass  densities  of 
©2  from  the  B-edge  to  the  entry  and  exit  points  of  the  target  cell. 

The  energy  available  for  hydro  due  to  Subgroup- 1 absorption  is 

DAIXU(l)  = PREFF(1,2)  X DAIX(l)  , 


where  PREFF(1,2)  is  the  pressure  efficiency  of  (Species  2)  absorbing 
Subgroup- 1 photons. 


For  the  target  cell,  the  effective  photon  fluence  for  Subgroup  1, 
consistent  with  the  energy  deposition,  is 


vTiTYm-  P DAIX(l)  X 2 X PMOXY 
ILUX(l)  ERGU(l)  X SIGUd,  2) 

^2 


p DAIX(l) 

|Q^|  ERGU(l)  X SIGUd, 2)  ’ 
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since 


1 ^ 2 X PMOXY 

[O2J  ^02 

For  Subgroups  2 through  5,  the  Subgroup- L energ-^  absorbed 


per  gram  is 


DAIX(L)  = - ERGU(L)  x SNR2(L)  ( -SUMIN(L)  _ ^-SUMOUT(L)J 

/ _ /o \ „3  \ I ) 


(p/3)  (R"  , - R"  ) 
out  in 


where 

4 

SUMIN(L)  = SPECIN(J)  X SIGU(L,  J) 

J=1 

4 

SUMOUT(L)  = 2 SPECOUT(J)  X SIGU(L,  J) 

J-1 

and  SPECIN(J)  and  SPECOUT(J)  are  the  areal  densities  of  the  J-species 
from  the  B-edge  to  the  cell  entry-point  and  exit-point. 

For  Subgroup  L,  the  average  pressure  efficiency  for  a target 
cell  is  defined  as 


PEFBAR(L)  - 


PREFF(L,  J)  X SPECIE(J)  x SIGU(L,  J) 
J=1 

4 

SPECIE(J)  X SIGU(L,  J) 

J=1 


The  energy  available  for  hydro  due  to  Subgroup  L is 
DAIXU(L)  = DAIX(L)  X PEFBAR(L) 
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For  the  target  cell  the  effective  fluence  for  Subgroup- L photons, 
consistent  with  the  energy  deposition,  is 


FLUX(L)  ERGU(L)  X XMU(L) 

J-1 


where 


XMU(L)  = 


SIGU(L,  J)  X SPECIE(J) 
J=1 


4 

j=i 


The  total  Group-U  energy  available  for  hydro  is  the  sum  of  that 
for  Subgroup- 1 and  that  for  Subgroups  2 to  5,  i.  e.  , 


DAIXU  = 2]  DAIXU(L) 
L=1 


Case  2 (B-edge  radius  (SRR(2))  greater  than  target-cell  exit- radius 
(SRN)) 


Inside  the  B-edge  the  internal  energy  is  set  by  the  expression 
DAIXU  = SPINT 


where  SPINT  is  set  in  Block  Data  BLOCKH.  The  HARC  code  used  a 


.12 


value  of  2.  07  x 10  erg/g,  which  we  have  reduced  by  a factor  of  10, 


.12 


since  2.07  x 10  erg/g  corresponds  to 
'g/eV 


— ^ ^ 2.  325X  10  g nitrogen  atom  = 30.  08  — — — - 

I.OX  lO-'^eriL/eV  mtro|!en  at 


om 


15.04 


eV 


charged  particle 
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where  [charged  particle]  = ■|[N'^]  + [e]|.  This  energy  per  charged  parti- 
cle corresponds  to  a temperature  of  about  10  eV  since  the  kinetic  energy 
per  particle  is  (3/2)  kT.  Thus,  we  have  assumed  that  the  particles 
inside  the  B-edge  have  a temperature  equal  to  nearly  1 eV.  We  are 
neglecting  the  difference  between  the  oxygen  and  nitrogen  mass,  but  this 
approximation  is  satisfactory  since  the  assumed  value  of  about  1 eV  is 
somewhat  arbitrary. 

We  have  no  need  for  the  Group-U  fluences  at  a point  inside  the 
B-edge  since  it  is  assumed  that  only  atomic  ions  exist  there. 


Case  3 (B-edge  radius  (SRR(2))  between  target-cell  entry-radius  (SR) 
and  exit- radius  (SRN)) 


To  obtain  the  specific  energy  available  for  hydro,  we  use  a 
volume- weighted  average  of  the  values  for  the  two  portions  of  the  cell 
inside  and  outside  the  B-edge.  For  the  cell  portion  inside  the  B-edge, 
the  specific  energy  is  equal  to  SPINT,  as  in  Case  2.  For  the  cell  poiiion 
outside  the  B-edge,  the  specific  energy  is  determined  in  a manner  simi- 
lar to  that  for  Case  1. 


Again,  we  treat  dissociating  and  ionizing  subgroups  separately. 


For  Subgroup  1 in  the  cell  portion  outside  the  B-edge,  we  have 


DAlX(l)  = 


WU(1) 

I 


(P'3)(R 


out 


- ^ ) 
edge 


i,  -XMU1x020Ut| 
,1  - e 


where  the  subscript  'edge'  denotes  B-edge  and  the  subscript  'out'  has 
the  same  meaning  as  in  Case  1. 

For  Subgroups  2 through  5 in  the  cell  portion  outside  the  B-edge, 
the  Subgroup- L energy  absorbed  per  gram  is 


1 

i 

i 

i 
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DAIX(L)  = - ERGU(L)  x SNR2(L)  _ ^-SUMOUT(L)l 

(p/3)  (R^  , - ) I ^ ( 

out  edge 

where  SUMOUT(L)  is  computed  just  as  in  Case  1. 

The  energy  available  for  hydro  due  to  Subgroup  L is  computed 
just  as  in  Case  1. 

The  total  Group-U  energy  available  for  hydro  in  the  cell  por- 
tion outside  the  B-edge  is,  as  in  Case  1,  the  sum  of  that  due  to  the  five 
subgroups. 

The  weighted  average  of  the  values  of  the  deposited  energy 
available  for  hydro  in  the  two  portions  of  the  cell  is  taken  to  be 


where 


W1  - |r^ 


- 0/1 


I 2 3 j 

R , - r:^ 

I out  in) 


For  the  cell  portion  outside  the  B-edge,  the  effective  photon 
fluences  for  Subgroup  1 and  for  Subgroups  2 through  5,  consistent  with 
the  energy  deposition,  are  computed  just  as  in  Case  1. 

9.4.2  Pressure  To  Be  Used  for  Hydro 

In  Subroutine  PHEAT,  the  pressure  computed  from  the  equation 
of  state  [p  = (y-l)pe]  by  using  the  total  specific  energy  available  for 
hydro  from  the  deposition  of  the  Group-X,  -H,  and  -U  energies,  is  cur- 
rently interpreted  as  the  total  pressure.  However,  an  electron  pressure, 
Pg  = implied  by  the  computation  in  Subroutine  PCHEM  where 

a common  temperature  is  found  for  the  electron  kinetic  temperature, 

1 

the  Ng  vibrational  temperature,  and  the  0(  D)-to-0(  P)  population  ratio. 
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It  is  possible  that  a negative  value  for  the  heavy-particle  temperature 
obtains  when  the  electron  pressure  is  subtracted  from  the  total  pressure. 
Our  interim  procedure  in  Subroutine  PCHEM  by  which  we  insure  no 
negative  value  for  the  heavy-particle  temperature  needs  to  be  reviewed 
and  improved. 

As  an  interim  measure  we  have  adopted  the  HARC  procedure 
for  limiting  the  amount  of  specific  energy  when  the  various  contributions 
are  added  to  obtain  SPINTN;  however,  the  limiting  term  has  not  been 
derived  and  needs  to  be  reviewed  and  improved.  Another  possible 
shortcoming  is  that  the  current  procedure  permits  adding  energy  to  a 
region  within  a B-edge  for  a previous  event. 

9.4.3  Flow  Chart  for  Subroutine  PHEAT 

Figure  21  is  a detailed  flow  chart  for  Subroutine  PHEAT. 
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SHN3  = (.SILN'P  5=  ^'^n  (,SR)3  $ GSIdi'  .-ir)  ^ nUKLOC  , 2)  > (SRN3  - SR3)'3 

DAlXXRrK.  r)  = li:,X  (lr  GSR)|  (T15XIN  - J I3XOUT) 


ILLIAG,  miliahz,ed  in  PL.IME  and  set  in  HKDGK,  is  0 il  R-edt^e  radius  SRH(2)  < SKN 
and  is  1 if  SRU(2)  SUN. 


lUFlAG  - 

r 
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SRR3  =.  [,SRR(2))^ 


Croup-U  Lnorgy  Deposition 


Coml»inc  line  intettral.s  from  PINT  and  LEDGE.  Compute  attenuation  factors  for 
Sub{;ruup  L enter*n*i  and  exiting'  cell. 


020UT  - 02INT  - 02T  $ 02rS  = 02(JUT  - lK)2rNT 
SPECOUT(l)  SN2IN'T  - SN2U  S SPECOUT(2)  = S021NT  - S02U 
SPECOUT(3)  SNINl  -SNT'  $ vSPKCOUTl4 ) = SOINT  - SOU 
SPlCTN'(l)  --  SPl-COUj-iU  - D-SN21  $ SPi:ClN(2)  = SPECOUT(2)  - DS02I 
SPECLN(3)  --  Si*ECOrT(3)  - pSNi  $ SPEC1N’(4)  = SPECOUT(4)  - DSOJ 

SUMIN(L)  = ^ SPrXINtJj  ^ SlGUtL..T)  $ SUMOUT(L)  = ^ SPECOUTU)  > SIGUfL,  J) 
J-1  J 1 


Compute  enero'  absorption  DAIXteri:-  n)  and  effective  photon  f!uence  P'LUX^photon/ 
rtr.^)  consistent  tin*  energy  depo.*^ition  for  each  Su>>group  I.  1 = 1,5.  Subgroup  1 is 
02'dis.sociatinR,  other  subcroups  arc  lomzing. 


G- 


T 


SRR(2)>  SR 

zjEz: 


SRR(2)  < SR 


PAlXd)  = [WL'(1)/C::rI  ((.■xp{-XMU!  • 02!\)  - exp(-X.\U'l  • 020UT)] 

UAIXd.)-  |[WUd.)  - l.RGUil.)  X SNR2d.)l  GSR; 'exp[-SL:MlN(  I.)]  - c.xp[-KL-MOUT(I.)]l 


( 


f, 

'0 


CO 


Fig.  21.  Flow  Chart  of  Subroutine  PHEAT. 
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Fig.  21.  (Continued). 
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10.  PROMPT  CHEMISTRY  ASSOCIATED  WITH  GROUP-X  AND 
GROUP-U  ENERGY  DEPOSITION  (SUBROUTINE  PCHEM) 


Subi’outine  PCHEM  first  combines  (a)  the  changes  in  species 
densities  due  to  debris  deposition,  computed  in  Subroutine  HPCHEM, 
with  (b)  the  previous  species  densities  of  a cell  and  then  computes  the 
ionization  and  species  densities  from  the  Group-X  and  Group-U  energy 
depositions. 

Subroutine  PCHEM  contains  numerous  compromises  and  ad- 
justments made  necessary  by  using  energy  deposition  models  that  are 
(a)  instantaneous  and  (b)  mismatched  with  the  late-time  chemistry 
modules. 

10.  1 DETERMINATION  OF  INITIAL  SPECIES  FOR  COMPUTING 

THE  PROMPT  CHEMISTRY  ASSOCIATED  WITH  GROUP-X 

AND  GROUP-U  ENERGY  DEPOSITIONS 

Before  computing  the  prompt  chemistry  associated  with  the 
energy  depositions  of  Groups  X and  U,  we  must  combine  the  changes  in 
species  densities  due  to  debris  deposition  with  the  previous  species 
densities. 

The  changes  in  the  species  densities  and  electronic  thermal 
and  e.xcitation  energy  are  stored  in  the  BUF2  array  bv  Subroutine 
HPCHEM.  These  variables  are  defined  in  Table  3. 

The  incorporation  of  these  changes  may  seem  trivial,  but  it  is 
complicated  by  the  facts  that  (a)  the  late-time  chemistry  for  the  HAG 
treats  all  molecular  ions  as  if  thev  are  NO^  and  (b)  it  is  necessary  to 

i i 

appropriately  convert  the  N2  and  produced  by  the  heavy-particle  de- 
position into  NO*.  This  conversion  is  made  so  as  to  conserve  nuclei 
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Table  3.  Inputs  to  Subroutine  PCHEM  from  BUF2  Array  for 
Those  Cell  Quantities  Updated  by  Subroutine 
HPCHEM  as  a Result  of  Debris  Deposition. 


I BUF2(I) 


1 

Change  in  [e], 

-3 

cm 

2 

Change  in  [N("^S)], 

-3 

cm 

3 

Change  in  [N(^D)j, 

-3 

cm 

4 

Change  in  [O], 

-3 

cm 

5 

Change  in  [N2], 

-3 

cm 

6 

Change  in  [O2], 

-3 

cm 

7 

Change  in  [N^  ], 

-3 

cm 

8 

Change  in  [O^], 

-3 

cm 

9 

Change  in  O excitation  energy  and  electron  thermal 
energy  per  new  electron  (eV/electron) 

10 

Not  used 

11 

Change  in  [N^], 

-3 

cm 

12 

Change  in  [O^], 

-3 

cm 

12 


of  a given  mass  number  and  electronic  state  but  with  only  partial  regard 
to  energy  considerations.  The  appropriate  mass  conversion  coefficients 
are  based  on  the  same  dissociative  recombination  products  of  N^,  O^, 
and  NO  used  in  the  later  treatment  of  Group-X  and  -U  energy  deposi- 
tions. The  assumed  dissociative  recombination  products  are  as  follows 
■;  Volume  IP  ; 


N2.e 


N(^S)  + N(^D) 


O^.e 


NO  + e 


0.  25  N(^S)  + 0.75  N(^D)  - O 


To  illustrate  this  conversion  process,  we  note  that  if  we  have 
a NO"^  molecular-ion  density  equal  to  [N^],  the  number  of  N("^S)  atoms 
we  get  from  dissociative  recombination  of  NO^  is  only  one-fourth  as 

many  as  from  dissociative  recombination  of  N,l.  Thus,  to  the  initial 

4 ^ 

density  of  N(  S)  we  must  add  a number  equal  to  0.  75[N„].  Similarly,  if 

we  have  a NO  density  equal  to  [O^]?  the  number  of  N(^S)  atoms  we  get 

from  dissociative  recombination  of  NO^  is  0.  25[NO'^]  w'hereas  we  would 

4 

get  no  N(  S)  atoms  from  the  dissociative  recombination  of  [ol  ].  Thus, 

....  4 ^ 

to  the  initial  density  of  N(  S)  we  must  subtract  a number  equal  to 

0.25[02J. 

4 2 

To  prevent  the  density  of  N{  S),  N(  D),  or  O from  becoming 
negative  as  a result  of  converting  N2  and  to  NO^,  we  appropriately 
reduce  either  the  N2  or  (produced  by  the  heavy  particles)  and  cor- 
respondingly increase  N2  or  O2.  Of  course,  the  electron  density  must 
be  reduced  to  maintain  charge  conseiwation  when  we  reduce  the  density 
of  either  N^  or  O^. 
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More  details  of  the  steps  in  converting  N2  and  O2  to  NO  are 
given  on  the  first  page  of  the  flow  chart  for  Subroutine  PCHEM  (Fig.  22). 

We  assume  that  [CO2]  and  [He]  are  reduced  by  the  debris  de- 
position in  proportion  to  the  [N2]  reduction. 

10.  2 IONIZATION  AND  SPECIES  DENSITIES  FROM  GROUP-X 
DEPOSITION 

Computation  of  the  species  densities  from  Group-X  deposition 
is  done  very  simply. 

1.  Use  the  specific  energy  absorbed  from  Group-X  deposition 

(DAIXX(erg/g))  computed  in  Subroutine  PHEAT  and  the  mass  density 

3 3 

(RHO(g/cm  ))  to  get  the  corresponding  energy  density  (DELEX(eV/cm  )). 

2.  Assume  35-eV  of  energy  is  expended  in  forming  an  ion- 

pair. 

3.  Distribute  the  ions  in  proportion  to  the  neutral  species, 
but  allow  no  more  than  one  ion  pair  per  neutral  particle.  Such  a pre- 
scription is  perhaps  unduly  simplistic  in  view  of  the  considerable  amount 
of  information  available  concerning  x-ray  and  electron  energy  degi’ada- 
tion  in  the  atmosphere  (see,  e.g.,  MS-75,  SK-75,  TA-74):  however,  no 
reliable  prescription  is  yet  available  to  predict  the  species  for  an  arbi- 
trary initial  mix  of  species. 

4.  Assume  the  molecular  ions  undergo  instantaneous  dissoci- 
ative recombination  in  accordance  with  the  prescription  in  Section  10.  1. 

5.  Assume  that  [CO2]  and  [He]  are  reduced  in  proportion  to 
the  [N2]  reduction. 
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Fig.  22.  Flow  Chart  of  Subroutine  PCHEM. 
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Sum  O and  N nuclei  foi  later  use  in  mabs-conscrvation  check. 


SUMOS  - OATOMS  - 2 ' 02S  ♦ FNUS  ^ (;i>S  ♦ AN»Or.l’S 
Sl’MNS  - KN4SS  . KN2US  ♦ 2 x KN2S  ^ KNOS  KNPS  « AMOLPS 


XL: 


Heset  PNU)l,ri,  DMOI-P2.  and  HMOLP3  (to  an  arbitrarv  lari', o iiumljcr)  and  UAMOl.P  ((u  0)  to  i ruvide  1 
for  pu.sMbihtv  that  a Cas.o-2  cell  follows  either  a Case-1  or  a Case-3  cell-  I)LlinUi->nvS  of  3 casesi  same 
as  in  PJIEAT. 

Case  1;  fi-edi^e  short  of  coll  entry-point. 

Case  2.  13-edge  beyond  cell  exit-point. 

Case  3;  )3-edge  wiihin  cell. 

For  Case-2  cell,  bypass  x-ray  deposition;  otiierwise,  start  on  Cases  1 and  3. 


©s. 


Start  Croup-X  Deposition 


Use  Grou|)-X  absorbed  specific  energy  DAI>LX(erg/g)  from  THKAT  and  density  flilOfg  cm^)  to  ret 
Group-X  energy  density  I>KLKX(eV/cm^).  AlPXfip  cni^)  - lon-pur  density.  Distnliute  ions  in  pro- 
portion to  neutral  species;  S.M.NKUT  - sum  of  neutral  .species.  Limit  ion  density  to  1 per  neutral. 

DELEX  ^ DAtXX  > RHO.  1.6E-12  S AIRX  = UF.U-iX  35 

SMNF.UT  = EN2S  + 02S  . ENOS  t ENISS  ♦ EN2DS  * OATOMS  S AXDST  = Afl’X  SMKEUT 


AXDST  = 1 


EN2PX  = AXDST  x ES2S  - 1)  S 02PX  •=  .MaxlO.  AXDST  x 02.s  - 1) 

ENOPX  = Ma.xfO,  AXDST  x EN'OS  - 1)  S EN-IST’X  - MaxtO. AXDST  > KN-ISS  - I) 

EN2DPX  = Max(0,  AXDST  x EN2DS  - 1)  $ OTX  M.ixiO.  AXDST  > OATOMS  - 1) 


Group-X  final  species  after  dissocnative  recombination. 

EN2X  = Max(l,EN2S  - FN2PX)  S 02X  = Mrix(l,02S  - 02I>X) 

ENOX  = Maxd.ENOS  - EN'Ol’Xl 

EN4SX  = .MaxU.ENTSS  - ENTSl’X  . EN2PX  . 0.25  ENOPX) 

EN2DX  - Max(!,EN2DS  - EN2DPX  • EN2PX  . 0.75  ENOPXI  S 1 NX  - EN4SX  . EN2nx 
OATOMX  - .Max(l,OATOM.S  - Ol'X  . 2 . 02)'X  . ENOPX  $ l.NPX  - ENPS  . F.NNSl'X  . EN2DPX 
OPX  = OPS  . OPX  S ENEX  - ENPX  . OPX  ♦ AMOLPS 


Sum  O and  N nuclei  for  mass-coiiser-vation  check.  Heduct  CO2  and  He  in  prop.iition  to  Nn  decrease 


SUMOX  . OATO.MX  ♦ 2 x 02X  . ENOX  . OI>X  . AMOLPS 
SUMNX  ■ EN4SX  . EN2DX  . 2 > EN2X  . EN'OX  • ENPX  . AMOLPS 
EHCOX  (SUMOX  - SEMO.S)  SUMOS  S HU’NX  ■-  fSLMNX  - SUMNS)  SlTtlN-S 
CO20X  u CO20S  .<  EN2X  EN2.S  $ HEX  HI  ^ • EN2X  EN2S 


Fig.  22.  (Continued). 


! 

Start  Group- U Opposition  for  Cases  I and  3 


Compute  attenuation  factors  jFSlGUfJ)]  fr->:n  products  of  Gr«>up-U  fluences  and  riviss  sections.  Distin- 
guish between  ionizing  (I>'(C»2I)  and  total  (f>IGUi2)|  factors  for  Gy  Compute  prubability  PKOI3SL'(J)  of 
Species  J .surviving  Group-U  deposition;  J 1,2,  3,4  for  N2,  O2.  N,  O. 


5^ 

FSIGU(J)  FLUXU.)  ^ SIGU{L,  J),  J=^l,4 

L--2 

FSIG21  ^ FSIGU(2)  S FSIGU(2)  - FSIG2f  ♦ FLUX(I)  ^ SIGU(i,2) 

PROIlSi;(J)  exp|Max(-F^;iGU(J),-70)l,  .1  = 1,4 
PRBS21  - exp[Ma.x(-FSlG21,-7U)l 

t 

Group-U  initial  species  after  dissociative  recombination;  FS-73  treatment  does  not  include  N(^D)  or  NO 
as  absorbers;  should  revise  simple  treatment  here. 


EN2i;  ^ F:N2X  X PROaSUil)  S 02U  - 02X  • PROBSU(2J 
02PU  - 02.x  X (1  - PRBS21)  S FS3M1  - FS!GUt3)  - FSIGU(l) 


Fig,  22.  (Continued). 


Ill 


Check  conser\’:»lion  of  O and  N nuclei.  Lot  CO2  i«n'l  He  bo  clostroyod  in  propoition  to  N2  doj't rurlom 


SL'MOK  OA'IOMK  * 2 > 02t’  • KNOU  » OlHI  * A\TO!  PS  * r>AMOU 
SUMNF  FNISF  * LN2DF  * 2 • r.N'2Ii  ♦ LNOU  • I.NPV  • AMOI.PS 
FRCO  (SltMOS  - SUMOF),  SUMOS  $ FKCN  - (FUMNvS  - FUMN  I- 
CO20U  ^ CO20X  ^ PROFkSnU)  5 HFU  = HKX  • PUOH.SU(l) 


lUMOLP 
) SPMNS 


Conchido.s  cumbmed  troatmonl  of  Cases  1 and  3.  T'«  roniploic  ('aco  l,  p..  to  70  To  c unplclo  Case  3 

determine  roll  prc:porties  by  volumn-weiphtinp  Iwo  porfion.s  inside  and  outside  B-«  dj.r. 


©IS 


sm{(2)  1 SH 


-1, 


") 


Fig.  22.  (Continued). 
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Fi^.  22.  (Continued). 
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As  one  of  the  mass  conservation  checks  during  development, 
we  have  compared  the  sums  of  the  O and  N nuclei  after  the  computation 
of  the  species  resulting  from  Group-X  deposition  with  the  sums  before 
the  deposition. 

If  the  target  cell  is  a Case- 2 cell  as  defined  in  Section  9.  4 and  j 

again  below  in  Section  10.3,  no  computation  is  made  of  species  resulting  j 

from  Group-X  deposition. 

10.3  IONIZATION,  SPECIES  DENSITIES,  AND  RESULTANT 
PRESSURE  FROM  GROUP- U DEPOSITION 

In  computing  the  species  from  Group-U  deposition,  we  consider 
three  cases,  just  as  we  did  in  Subroutine  PHEAT: 

Case  1:  B-edge  radius  (SRR(2))  is  smaller  than  the  target-cell  entry 
radius  (SR). 

Case  2:  B-edge  radius  (SRR(2))  is  greater  than  the  target-cell  exit 
radius  (SRN). 

Case  3:  B-edge  radius  (SRR(2))  is  between  the  target-cell  entry  radius 
(SR)  and  the  exit  radius  (SRN). 

The  total  attenuation  factors  (FSIGU(J))  are  computed  from  the 

2 

summed  products  of  Group-U  fluences  (FLUX(L),  photons/cm  ) and 

2 

cross  sections  (SIGU(L,J),  cm  ).  (Recall  that  the  Group-U  fluences 
computed  in  Subroutine  PHEAT  are  effective  fluences,  consistent  with 
the  absorbed  energy  computed  as  the  difference  of  the  (steradianal)  en- 
ergy entering  and  leaving  the  cell.)  The  attenuation  factors  are 

5 

FSIGU(J)  = J]  FLUX(L)  X SIGU(L,  J)  , J = l,4 
L-2 

FS1G2I  = FSIGU(2) 
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FSIGU(2)  - FSIG2I  + FLUX(l)  X SIGU(1,2)  , 

where  SIGU(L,  J)  is  the  absorption  cross-section  of  Species  J for 
Subgroup- L photons,  tabulated  in  Table  4 and  entered  as  data  in  Block 
Data  BLOCKH,  and  where  we  have  distinguished  between  the  ionizing 
(FSIG21)  and  total  (FS1GU(2))  attenuation  factors  for  02* 

The  total  survival  probability  of  Species  J is 

PROBSU(J)  = . 

The  probability  of  not  being  ionized  is 

PRBS21  = e-FSIG21 


If  the  species  densities  after  the  Group-X  deposition  (but  before 

the  Group-U  deposition)  are  denoted  by  a subscript  x,  the  initial  species 

after  the  Group-U  deposition,  with  provision  for  instantaneous  dissoci- 

2 

ative  recombination  of  the  molecular  ions  but  without  production  of  N(  D), 
are  [FS-73J 


[Njlu  = PROBSU(I) 

i°2*u  " PROBSU12) 


lN|^=  [N|,xi>ROBSU(3).2x1N2^x  fSIGU(1)x£^W[^-|^^ 
101„  = lO^^x  PR0BSU(4)  . 2X  lOj^x  FSIGU(2)  x TsiGUH)  I FsS'* 

[N\=  In'|^-2{1N2J^-|N2U.  [Nl^-lNl^ 
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Table  4.  UV  Absorption  Cross-Sections  SIGU(L,  J)^ 
(in  units  of  10”1^  cm^). 


Species 

J 

SIGU(L,  J) 

L = 1 
11.0  eV 

L = 2 
14.  8 eV 

L = 3 
16.  2 eV 

L = 4 
22.0  eV 

L = 5 
36.  0 eV 

^2 

1 

0 

8 

20 

25 

20 

2 

2 

16 

20 

25 

20 

N 

3 

0 

10 

10 

10 

10 

O 

4 

0 

3.  2 

3.  2 

9 

9 

^From  Table  A-1  of  FS-73,  in  turn  from  DD-70b. 


Since  the  Group-U  deposition  procedure  does  not  provide  for 
2 

either  N(  D)  or  NO,  as  an  interim  measure  we  partition  the  N atoms  be- 
4 2 

tween  N(  S)  and  N(  D)  according  to  the  prescription 
[N(^S)1^  = [Nl^  [N(^S)1^/[N1^ 

- [Nl^]l  - LN('^S)l^/[Ny  , 

decrease  the  NO  in  proportion  to  the  N^  decrease,  and  increase  N and 
O correspondingly: 

[NOl^  - [NO]^  X PROBSU(l) 

6[N0]  = [NOl^{l  - PROBSU(l)} 

[N(‘^S)1^  - [Nl'^S)]^  + 6[N0] 

[Nl^  -»  [Nl^  . 6[NOl 
[Ol^  - [Ol^  4 6[NO]  . 
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To  maximize  the  ionization  at  early  time,  as  suggested  by 
W.  Knapp,  we  re-establish  the  molecular  ions  produced  and  dissoci- 
atively  recombined  by  the  Group-X  and  -U  depositions  (and  allow  the 
molecular  ion  decay  to  be  predicted  by  the  later-time  chemistry  module). 
The  increment  in  such  molecular  ions  is 

DAMOLP  s 6M^  = [Njl^  + [0^|^  + iNO^l^  + 

* 

where 

The  re-establishment  of  all  molecular  ions  as  having  the  properties  of 

4-  -f- 

NO  requires  subtraction  of  all  the  products  of  NO  dissociative  recom- 
bination from  the  appropriate  species  so  as  to  conserve  N and  O nuclei: 

[N(^S)]^  - [N(^S)]^  - 0.  25  X 6m^ 

[N(^D)]^  = [Nt^D)|^  - 0.  25  X 6M^ 


|Olj  = [O]^-  6m"  . 

However,  if  [N("^S)lj,  |N(^D)]j,  or  [Oj^  is  negative,  we  must  reduce 
so  that  none  of  them  is  negative. 


The  final  values  for  the  species  densities  after  Group-X  and  -U 
depositions  and  re-establishment  of  the  molecular  ions  are  1^2 1^- 
[NOl^,  |N(‘^S)1^,  [N(^D)|^,  [Olj,  [n'1^,  lO^j^,  and  [e].  the  electron 
density, 

|e|  = [N^l^  ^ [0"1^  . 6m‘  . m; 
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where  is  the  molecular-ion  density  prior  to  Group-X  deposition. 

There  remains  the  task  of  devising  a truly  satisfactory  pro- 
cedure for  computing  a pressure  increment  in  addition  to  that  computed 
in  Subroutine  PHEAT.  The  treatment  in  Subroutine  PHEAT  is  incom- 
plete because  the  Group-U  subgroups  are  deposited  independently  of  each 
other  across  an  appreciable  path  length  without  accounting  for  either 
species  depletion  by  other  subgroups  or  additional  absorbers  resulting 
from  the  (instantaneous)  dissociative  recombination  that  is  included  in 
Subroutine  PCHEM.  (Note,  however,  that  the  partial  heating  computed 
in  Subroutine  PHEAT  does  serve  the  purpose  of  defining  effective  flu- 
ences  for  the  Group-U  subgroups  that  are  used  in  Subroutine  PCHEM.) 
Any  additional  pressure  increment  computed  in  Subroutine  PCHEM  must 
be  consistent  with  (1)  that  computed  in  Subroutine  PHEAT  and  with 
i (2)  the  return  of  pressure  computed  by  the  later-time  chemistry  in 

\ CHEMEF.  The  problem  is  further  complicated  because  the  molecular 

[ ions  which  are  produced  by  the  Group-X  and  -U  depositions  and  which 

undergo  dissociative  recombination  as  part  of  the  modeling  in  Subroutine 

PCHEM  are  N,!,  and  ol  whereas  the  molecular  ions  that  are  re- 
[ 2 2 

: established  are  NO*^  ions  which  are  carried  in  CHEMEF.  As  an  interim 

measure  we  set  the  pressure  increment  DEEP  equal  to  zero. 

I To  ensure  that  the  heavy-particle  pressure  will  always  be 

: positive,  we  have  added  the  electron  pressure  (computed  below)  to  the 

previously-found  value  for  the  total  pressure.  This  interim  measure 
needs  to  be  reviewed  and  improved. 

To  compute  the  temperature  TEE  taken  to  be  common  to  tiie 

1 3 

thermal  electrons,  the  N2  vibrational  states,  and  the  0(  D)-to-0(  P) 
population  ratio,  we  use  the  GET-prepared  Subroutine  TEXK  Vol.  IT 
which  determines  such  a temperature  when  given  " O’’ , 
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[e],  and  the  energy  per  cubic  centimeter  (ET)  to  be  divided  among  the 
three  modes.  As  an  interim  measure  we  compute  ET  from 

ET  = 2 X 6[el  f 6[e  L„  x E ^ 

where  ^ is  the  incremental  electron  density  produced  by  the 

Group-X  and  -U  depositions,  the  factor  2 is  an  interim  assumption  for 

the  energy  in  the  three  modes  per  newly-formed  electron  in  the  Group-X 

and  -U  depositions,  is  the  incremental  electron  density  produced 

by  the  heavy-particle  depositions,  and  E is  the  sum  of  the  electron 

pe 

thermal  and  0(  D)  excitation  energies  per  newly-formed  electron  in  the 
heavy-particle  deposition,  obtained  from  Subroutine  HPCHEM. 

As  another  of  the  mass  conservation  checks  during  develop- 
ment, we  have  compared  the  sums  of  the  O and  N nuclei  after  the  com- 
putation of  the  species  resulting  from  Group-U  deposition  with  the  sums 
before  the  deposition. 

As  with  the  Group-X  deposition,  we  assume  as  an  interim 
measure  that  ICO^]  and  [He]  are  reduced  in  proportion  to  the  [N2] 
decrease. 

The  foregoing  describes  the  combined  treatment  of  target  cells 
that  are  Case  1 or  Case  3.  For  Case  1,  the  treatment  is  complete,  but 
for  Case  3,  the  final  cell  properties  are  determined  by  volume-weighting 
the  two  portions,  one  of  which  is  inside  the  B-edge  and  the  other  is  out- 
side the  B-edge. 

For  Case  2,  the  entire  cell  is  inside  the  B-edge  and  the  only 
remaining  species  are  assumed  to  be  N , O^,  and  e.  The  temperature 
of  both  electrons  and  heavy  particles  is  taken  to  be  1 eV. 
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The  final  step  in  Subroutine  PCHEM  is  to  store  the  cell  quanti- 
ties in  temporary  storage  BUFl  until  the  loop  over  the  entire  column  of 
target  cells  is  completed,  after  which  the  entire  column  of  cell  quantities 
is  transferred  to  scratch  storage  in  large-core  memory. 
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